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ABSTRACT 


An  experimental  wind  tunnel  investigation  has  been  conducted  to 
determine  the  hydrodynamic  characteristics  of  bluff  symmetrical 
sections  with  high  thickness/chord  ratios.  The  sections  are  used  as 
fairings  for  circular  cylindrical  members,  such  as  towing  cables  and 
offshore  drilling  rig  riser  pipes,  which  are  deeply  immersed  in  the 
ocean  environment  so  that  they  are  cavitation  free.  The  fairings 
serve  to  reduce  drag  and  lateral  vibrations  due  to  vortex  shedding. 
The  section  tested  had  a  40  percent  thickness/chord  ratios. 
Measurements  included  the  two-dimensional  minin.um  drag  coefficient, 
chordwise  neutral  stability  point,  yaw  torque  about  the  pivot  center, 
and  boundary  layer  transition  and  separation  location. 

This  research  was  carried  out  under  the  Naval  Sea  Systems  Command 
General  Hydromechanics  Research  Program  administered  by  the  David  W. 
Taylor  Naval  Ship  Research  and  Development  Center,  under  Office  of 
Naval  Research  Contract  No.  N0014-82-K-0006. 

1.0.  INTRODUCTION 

Marine  applications  of  line  structures  (high  length/diameter 
ratio)  with  circular  cylindrical  sections  include  moored  and  towed 
cable  systems,  risers  for  offshore  drilling  rigs,  and  the  "cold  water 
pipe"  for  the  Ocean  Thermal  Energy  Conversion  (OTEC)  system.  All  of 
these  systems  experience  relative  motion  between  the  line  structure 
and  the  surrounding  water,  due  either  to  ship  motion  or  current  and 
wave  motion.  Because  of  this  relative  motion,  the  line  structure 
experiences  unsteady  (time-dependent)  hydrodynamic  forces  due  to 
vortex  shedding.  When  the  frequency  of  the  vortex  shedding  is  close 
to  the  natural  frequency  of  the  line  structure,  a  resonant  structural 


response  condition  will  occur.  This  resonant  condition  will  increase 
the  drag  force  (force  parallel  to  current  field)  which  in  turn  will 
result  in  cyclic  motions  and  stresses. 

One  solution  to  the  shed  vortex  problem  is  to  fair  the  circular 
cylinder  with  a  streamlined  shape.  An  obvious  drawback  is  that  the 
fairing,  unless  it  is  allowed  to  swivel  freely  and  align  itself  with 
the  flow  field,  will  act  as  a  wing  and  generate  large  transverse  lift 
forces.  However,  this  in  turn  dictates  that  the  position  of  the 
hydrodynamic  (aerodynamic)  center  must  be  aft  of  the  mechanical  center 
of  rotation  (central  axis  of  the  line  structure)  for  weathervane 
stability. 

It  should  be  noted  that  the  hydrodynamic  center,  or  neutral 
stability  point,  is  used  in  contrast  to  the  center  of  pressure.  The 
hydrodynamic  center  is  defined  as  the  position  along  the  chord  about 
which  the  moment  coefficient  is  constant  independent  of  lift 
coefficient,  while  the  center  of  pressure  position  varies  with  lift 
coefficient.  For  symmetrical  sections,  this  constant  moment 
coefficient  has  a  value  of  zero.  The  neutral  point  is  defined  as  the 
point  along  the  chord  where  the  slope  of  pitching  moment  coefficient 
versus  lift  coefficient  curve  is  zero.  Therefore  the  terms 
hydrodynamic  center  and  neutral  point  are  synonymous  for  a  symmetrical 
section  and  are  used  interchangeably. 
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2.0.  PROBLEM  DESCRIPTION 

Rigid  symmetrical  fairings  in  short  lengths  are  presently 
commercially  available  from  Fathom  Oceanology  Ltd.  for  application  to 
towing  cables,  Fig.  2.1.  The  fairings  range  in  size  from  chord 
lengths  of  5  to  15  cm  for  towing  cable  applications,  to  2  m  for  a 
drill  pipe  riser  fairing,  [1].  Hydrodynamic  problems  with  the  towing 
cable  fairings  have  been  encountered,  as  discussed  by  Henderson  [2]. 
A  fairing  section  which  had  a  thickness/chord  ratio  (t/c)  of  25 
percent  was  found  to  have  a  hydrodynamic  center  position  at  the 
15.4  percent  chord  position,  which  was  aft  of  the  center  of  rotation 
by  a  distance  of  only  2.3  mm  (2.9  percent  chord).  It  was  determined 
experimentally  that  boundary  layer  separation  over  the  aft  portion  of 
the  section  caused  this  problem.  The  addition  of  a  flat  plate 
trailing  edge  extension,  to  fix  the  position  of  the  aft  stagnation 
point,  increased  the  chord  to  100  mm  and  resulted  in  the  movement  of 
the  hydrodynamic  center  to  the  25  percent  chord  position.  The  large 
diameters  of  line  structures  such  as  drill  riser  pipes  (about  6.5  to 
10  m)  require  that  the  fairing  chord  length  be  as  small  as  possible 
for  handling  and  installation  considerations.  This  requirement 
results  in  fairing  sections  with  very  high  thickness  to  chord  radios 
(up  to  50  percent).  Boundary  layer  separation  was  apparently  also 
observed  in  a  series  of  wind  tunnel  tests,  Grant  [3],  on  high 
thickness/chord  ratio  (0.4  to  0.5  c)  fairings  designed  for  riser 
lines.  It  was  observed  that  the  fairing  did  not  possess  weathervane 
stability.  In  this  particular  application,  the  problem  was  solved  by 
the  addition  of  fins  to  the  fairing  trailing  edge,  which  acted  like  a 
split  flap  and  stabilized  the  fairing. 
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3.0.  PROBLEM  APPROACH 


The  selection  of  an  optimum  airfoil  section  shape  for  use  as  a 
fairing  must  include  consideration  of  the  following: 

(1)  streamlined  symmetrical  section  for  low  drag 

(2)  position  of  maximum  thickness/chord  as  close  to  the 

leading  edge  as  possible  in  order  for  center  of 

rotation  to  be  forward  of  hydrodynamic  center  for 
weathervane  stability 

(3)  high  thickness/chord  ratio  (bluff)  sections  to  reduce 
the  chord  size  of  the  fairings 

(4)  separation  free  boundary  layer. 

The  utilization  of  bluff  sections  as  fairings  will  lead  to 
problems  with  boundary  layer  separation,  as  already  noted.  Some  means 
must  therefore  be  used  to  prevent  or  control  boundary  layer  separation 
and  to  provide  rotational  stability.  Three  techniques  present 

themselves  as  appropriate  solutions,  either  separately  or 
synergistically .  These  include  the  use  of  separation  resistant 
sections,  vortex  generators  for  boundary  layer  control  and  trailing 
edge  wedges  for  stability. 

Recent  advances  in  the  development  of  boundary  layer  separation 
resistant  airfoil  sections  with  specified  pressure  profiles  have  shown 
some  interesting  results.  A  section  with  a  53.6  percent  t/c,  designed 
by  Liebeck,  Smith  [4],  is  the  result  of  shaping  the  forward  position 
of  the  section  so  that  the  pressure  gradient  is  favorable  to  laminar 

flow,  and  then  using  the  Stratford  pressure  recovery  over  the  rear 

portion  for  separation  control. 

A  passive  technique  for  boundary  layer  control  is  the  use  of 
vortex  generators.  This  technique  relies  on  the  increased  mixing 


between  the  external  stream  and  the  boundary  layer  as  promoted  by 
vortices  trailing  longitudinally  over  the  surface  from  the  generators. 
Vane  type  generators  are  the  ones  most  often  used.  They  consist  of  a 
row  of  small  plates  that  project  normal  from  the  surface  with  each  one 
set  at  an  angle  of  incidence  to  the  local  flow  to  produce  a  single 
trailing  vortex. 

The  use  of  "split-flap"  wedges  on  the  trailing  edge  of  rudders 
has  been  investigated,  Thieme  [5],  as  a  means  to  control  the 
hydrodynamic  moment.  This  suggests  itself  as  a  technique  which  might 
have  application  for  the  fairing  configuration  to  achieve  rotational 


4.0.  SYMMETRICAL  FAIRING  SECTIONS 
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4.1.  Section  Geometry 

A  total  of  five  symmetrical  fairing  sections  were  selected  for 
this  study.  These  included: 

(1)  Fathom  fairing  section 
NACA  0040 
Liebeck 
JfS  61-TR-40 
JfS  62-TR-40. 

A  thickness/chord  ratio  (t/c)  of  40  percent  was  arbitrarily  chosen  as 
the  maximum  that  might  be  successful  in  minimizing  the  chord  length 
while  hopefully  operating  with  a  separation  free  boundary  layer.  The 
one  exception  was  the  Fathom  fairing,  with  a  t/c  of  25  percent.  This 
section  was  selected  to  provide  a  baseline  for  data  comparisons  with 
past  studies  of  this  section,  and  to  measure  the  degree  of 
improvement,  if  any,  provided  by  the  thicker  sections. 


(2) 

(3) 

(4) 

(5) 


m 

The 

section  geometry  may  be  described  primarily  by 

.y 
’  */ 

(Fig.  4. 

1)  . 

;y 

(1) 

maximum  thickness  (t/c) 

(2) 

chordwise  position  of  maximum  thickness  (xt/c) 

A 

(3) 

nose  radius  (r/c) 

(4) 

trailing  edge  thickness  (tg/c) . 

In  addition,  the  shape  of  the  section  aft  of  the  maximum  thickness 
location  may  be  described  as  either  concave  or  convex,  Fig.  4.1.  The 
individual  fairing  geometry  is  contained  in  Table  4.1,  the  sections 
are  shown  in  Fig.  4.2,  and  the  coordinates  contained  in  Table  4.2. 
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Fig.  4.1.  Definition  of  geometric  variables. 
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Table  4.1  -  Pairing  Geometry 


Section 

t/c 

Xt/C 

r/c 

After  Section 

Fathom 

0.25 

0.125 

0.125 

Convex 

NACA  0040 

0.40 

0.30 

0.125 

Convex 

Liebeck 

0.40 

0.35 

0.1763 

Concave 

JfS  61-TR-40 

0.40 

0.20 

0.25 

Concave 

JfS  62-TR-40 

0.40 

0.20 

0.25 

Concave 

4.2.  Section  Description 
4.2.1.  Fathom  Flexnose  Fairing 

The  patented  Fathom  Flexnose  fairing  is  one  of  the  few 
commercially  available  cable  fairings  and  is  sold  by  Fathom  Oceanology 
Ltd.  of  Port  Credit,  Ontario,  Canada.  The  fairing  is  available  in 
t/c's  ranging  from  0.13  to  0.27  and  chord  lengths  of  6.35  cm  to 
42.44  cm.  Fig.  2.1  shows  fairings  of  chord  lengths  of  6.35  and 
12.19  cm. 

Henderson  [2]  has  determined  the  hydrodynamic  and  mechanical 
characteristics  of  the  Fathom  Flexnose  25  percent  fairing.  Lift, 
drag,  and  pitching  moment  coefficients  were  measured  at  a  Reynolds 
number  of  2.2  x  10s ,  which  corresponds  to  towing  a  5.08  cm  fairing  at 
9.8  knots.  Qualitative  boundary  layer  flow  visualization  studies  were 
also  done  at  this  Reynolds  number.  The  aerodynamic  center  was 
determined  from  the  lift  and  pitching  moment  data.  A  modification  to 
the  Fathom  fairing  trailing  edge  was  made  in  an  attempt  to  improve  its 
hydrodynamic  behavior.  The  two-dimensional  hydrodynamic  drag  polars 
were  also  measured  for  this  modified  fairing. 

Fathom  Oceanology  Ltd.  produces  several  Flexnose  fairings  of 
slightly  varying  profiles.  The  profile  chosen  for  this  investigation 


Lie beck  40% 

NACA 

Fathom 

JfS  62M 

JfS  61M-TR-40 

0040 

Flexnose 

TR-40 

x/c 

+y/c 

x/c 

+y/c 

+y/c 

+y/c 

x/c 

+y/c 

0.00000 

0.00000 

0.000 

0.00000 

0.00000 

0.00000 

0.000 

0.00000 

0.01075 

0.04319 

0.004 

0.03654 

0.03137 

0.03980 

0.004 

0.03980 

0.03252 

0.08273 

0.010 

0.05679 

0.04899 

0.06245 

0.010 

0.06245 

0 . 06626 

0.11841 

0.020 

0.07866 

0.06782 

0.08718 

0.020 

0.07866 

0.10313 

0.14402 

0.040 

0.10759 

0.09165 

0.12000 

0.040 

0.12000 

0.15717 

0.16874 

0.060 

0.12792 

0.10677 

0.14283 

0.060 

0.14283 

0.22017 

0.18653 

0.080 

0.14357 

0.11662 

0.16000 

0.080 

0.16000 

0.27713 

0.19561 

0.120 

0.16627 

0.12490 

0.18330 

0.120 

0.18330 

0.33586 

0.19874 

0.125 

0.16851 

0.12500 

0.18540 

0.125 

0.18540 

0.40512 

0.19310 

0.160 

0.18139 

0.12500 

0.19596 

0.160 

0.19596 

0.44298 

0.18154 

0.200 

0.19125 

0.12494 

0.20000 

0.200 

0.20000 

0.49094 

0.15644 

0.250 

0.19804 

0.12470 

0.270 

0.19467 

0.56479 

0.11647 

0.275 

0.19958 

0.12448 

0.19381 

0.340 

0.17957 

0.65158 

0.07907 

0.300 

0.20006 

0.12417 

0.410 

0.15709 

0.74160 

0.04832 

0.350 

0.19829 

0.12323 

0.17607 

0.480 

0.13040 

0.82670 

0.02564 

0.425 

0.19002 

0.12080 

0.14931 

0.550 

0.10299 

0.89994 

0.01099 

0.500 

0.17647 

0.11680 

0.11730 

0.620 

0.07824 

0.96434 

0.00233 

0.575 

0.15878 

0.11083 

0.08425 

0.690 

0.05891 

1.00000 

0.00056 

0.650 

0.13775 

0.10250 

0.05421 

0.760 

0.04648 

0.725 

0.11388 

0.09142 

0.03064 

0.830 

0.04093 

0.800 

0.08744 

0.07718 

0.01546 

0.090 

0.04000 

0.875 

0.05848 

0.05940 

0.00864 

1.000 

0.04000 

0.950 

0.02689 

0.03769 

0.00750 

1.000 

0.00420 

0.02083 

0.00750 

Table  4 . 

2.  Non-Dimensional 

Coordinates  of  the  Tested  Fairing 

Sections 

was  the  same  as  the  unmodified  profile  that  was  used  in  Henderson's 
study.  This  fairing  was  named  Flexnose  B  by  Wingham  [6]. 

The  experimental  Reynolds  number  range  was  extended  from 
1.0  x  106to  6.5  x  106to  complete  the  picture  of  the  behavior  of  the 

Flexnose  25  percent  fairing.  The  work  of  Henderson  could  thus  be  used 
as  a  basis  for  comparing  with  the  results  of  this  study. 

4.2.2.  NACA  0040  Section 

The  characteristics  of  the  NACA  00XX  series  of  symmetrical 
airfoils  are  well  documented.  Goett  and  Bullivant  [7]  tested  the  NACA 
0009,  0012,  and  0018  airfoils  in  1938.  In  1940,  Bullivant  [8] 
continued  the  study  with  the  NACA  0025  and  0035  sections.  Summaries 
of  the  results  of  the  NACA  00XX  tests  are  contained  in  Eastman,  Jacobs 
and  Abbott  [9],  Abbot  and  von  Doenhoff  [10],  Hoerner  [11],  and 
Carmichael  and  Meggitt  [12],  All  of  the  airfoil  characteristics  for 
these  sections  were  measured  at  high  Reynolds  numbers,  ranging  from 
106  to  107.  Althaus  [13]  tested  the  NACA  0033  section  at  Reynolds 
numbers  of  8.0  x  104,  1.2  x  105,  and  1.5  x  10s.  Although  these 

Reynolds  numbers  are  at  the  low  end  of  the  range  of  present  interest, 
the  reported  drag  coefficients  indicated  that  an  NACA  00XX  section 
might  be  an  acceptable  choice  for  a  cable  fairing. 

The  NACA  0040  section  was  thus  chosen  for  testing  in  part  because 
this  extensive  set  of  data  from  earlier  work  could  also  act  as  a  basis 
of  comparison  with  the  results  the  present  investigation.  In 
addition,  it  also  serves  as  an  example  of  a  section  with  a  convex 
after  section. 

4.2.3.  Liebeck  Section 

The  Liebeck  fairing  [4]  was  designed  such  that  laminar  flow  would 


be  maintained  along  the  forebody  with  no  boundary  layer  separation 
occurring  in  the  pressure  recovery  region  at  a  Reynolds  number  of  107. 
The  maximum  t/c  that  theoretically  met  those  requirements  was  53.6 
percent.  The  shape  of  tne  aft  section  was  designed  using  the  theory 
of  Stratford  [14]  for  an  imminent  separation  pressure  recovery 
velocity  distribution.  Thus  the  afterbody  was  termed  a  "Stratford 
pressure  recovery  section."  For  this  investigation  the  Liebeck  fairing 
was  reduced  to  a  40  percent  t/c  by  a  linear  reduction  of  the 
coordinates.  This  shape  was  chosen  for  the  study  primarily  because  of 
the  predicted  low  drag  for  its  high  t/c,  and  as  an  example  of  a 
section  with  a  convex  after  section. 

4.2.4.  JfS  Sections 

JfS  sections  [5]  are  members  of  a  family  of  balanced  ship  rudder 
profiles  developed  at  the  Shipbuilding  Institute  of  the  University  of 
Hamburg.  These  profiles  were  developed  by  a  systematic  variation  of 
form  parameters  which  included:  cross-sectional  area,  leading  edge 
radius  of  curvature,  length  of  forebody,  length  of  pressure  recovery 
section,  length  of  trailing  edge,  thickness  of  trailing  edge,  surface 
slopes,  and  radius  at  maximum  thickness.  The  profile  coordinates  are 
expressed  as  polynomials  composed  of  so-called  influence  functions. 
One  influence  function  was  written  for  each  of  the  form  parameters. 
Thieme  [5]  gives  a  complete  description  of  this  technique.  Two  JfS 
sections  of  thickness/chord  ratios  of  25  percent  were  tested  by 
Thieme,  the  JfS  61-TR-25,  and  the  JfS  62-TR-25.  Lift,  drag,  and 
pitching  moment  coefficients  were  measured  at  Reynolds  numbers  of 
0.56  x  10fc,  0.70  x  106,  and  0.79  x  106 . 

Two  modifications  were  made  to  the  JfS  sections  for  this  project. 
The  thickness  ratio  was  increased  from  20  to  40  percent,  and  the 
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leading  edge  radius  was  modified  to  a  forebody  which  was  circular  in 
section  to  the  maximum  thickness  position.  Because  of  these 
modifications,  the  JfS  sections  were  renamed  the  JfS  61M-TR-40  and 
JfS  62M-TR-40. 

4.3  Circular  Cable  Fairing  Efficiency 

An  important  geometrical  consideration  for  a  symmetrical  section 
designed  to  fair  a  circular  cylindrical  cable  is  that  the  circular 
diameter  be  inscribed  in  the  fairing  as  far  forward  as  possible. 
Fig.  4.3  shows  the  nondimensional  diameter,  d/c,  of  these  circles  as  a 
function  of  chord  location.  Each  circle's  center  corresponds  to  the 
mechanical  pivot  location  of  the  fairing.  The  largest  diameter  circle 
for  the  NACA  0040  is  40  percent  of  the  chord  located  at  x/c  =  0.3, 
39.6  percent  located  at  x/c  =  0.35  for  the  Liebeck  section  and  40.0 
percent  located  x/c  *0.2  for  the  JfS  sections.  Thus  the  JfS  sections 
are  the  most  efficient  in  terms  of  the  packaging  of  the  circular 
cylinder . 

4.4.  Fairing  Surface  Slope 

An  important  geometrical  characteristic  of  the  fairing  shape  is 

the  local  surface  angle.  Since  the  pressure  gradient  is  a  function  of 

this  angle,  the  surface  slopes  for  each  of  the  sections  as  a  function 

of  the  chord  are  shown  in  Fig.  4.4.  Also  indicated  are  the  trailing 

edge  surface  angle,  0  .  Note  that  both  the  NACA  and  Fathom  fairings 

te 

have  convex  pressure  recovery  sections,  while  the  Liebeck  and  JfS 
sections  have  concave  pressure  recovery  sections.  For  both  JfS 
sections,  the  maximum  slope  in  the  pressure  recovery  section  is  less 
than  the  maximum  slope  for  the  Liebeck  section. 
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5.0.  EXPERIMENTAL  EQUIPMENT 


5.1.  Wind  Tunnel 

The  use  of  a  wind  tunnel  for  characterizing  the  fairing  section 
hydrodynamics  is  appropriate,  as  long  as  Reynolds  number  equivalence 
is  maintained.  If  it  is  assumed  that  the  fairing  section  is  being 
operated  at  depths  great  enough  to  avoid  cavitation,  the  aerodynamic 
characteristics  will  be  the  same  as  the  hydrodynamic  characteristics. 
The  tests  were  conducted  in  the  University  of  Washington  "Venturi" 
open  return  wind  tunnel,  which  has  an  octagonal  cross-section  of  0.9  m 
across  vertices  and  0.79  m  across  flats.  With  a  maximum  speed  of 
27  m/s,  a  fairing  section  model  with  a  30.5  cm  chord  may  be  operated 
at  a  Reynolds  number  up  to  5  x  105.  The  operational  Reynolds  number 
range  for  faired  cables  is  from  5  x  104  to  1  x  106,  while  it  is  about 
2  to  5  x  10fa  for  faired  riser  pipes. 

5.2.  Tunnel  Turbulence  Intensity 

The  wind  tunnel  turbulence  intensity  was  measured  using  the 
sphere  drag  coefficient  method  outlined  in  Pope  [15],  "Wind  Tunnel 
Testing."  The  sphere  drag  coefficient  was  measured  as  a  function  of 
diameter  Reynolds  number.  The  Reynolds  number  at  which  the  drag 
coefficient  equals  0.30  is  termed  the  "critical  Reynolds  number,"  Rcrr 
which  has  a  value  of  385,000  in  free  air.  The  turbulence  factor  is 
then  defined  as 

TF  =  385,000./Rcr.  (5.1) 
For  the  Venturi  Wind  Tunnel,  Rcr  was  equal  to  1.5  xlO5,  which  is 
equivalent  to  a  TF  equal  to  2.85.  The  turbulence  intensity  level,  T^, 
is  defined  as 

Tu  =  u/Uro  (5.2) 
where  u  =  root  mean  square  (rms)  velocity  fluctuations 


Um  =  free  stream  mean  velocity. 

The  turbulence  intensity  is  related  to  the  turbulence  factor  so  that 
=  2.15  percent  for  the  Venturi  Wind  Tunnel. 

Since  atmospheric  turbulence  levels  in  free  flight  are  very  much 
less  than  one  percent,  it  is  common  to  take  means  to  reduce  Tu  to 
values  below  this.  However,  this  problem  must  be  reassed  for 
conditions  in  the  ocean.  In  general,  the  turbulence  intensity  will 
vary  with  location  and  depth.  Because  of  the  difficulties  of 
measuring  the  ocean  turbulence  level,  it  is  common  to  use  a  parameter, 
£,  which  is  the  rate  of  viscous  dissipation  of  velocity  fluctuations 
Hinze,  [16].  The  rms  velocity  fluctuations  are  then  defined  as 


/  [A(e)2N3(k.f  ^dk 
ki  1 


(5.3) 


where 


=  wave  number 
*  0.1  cycle/m 
=  5.0  cycle/m 


=  constant  =  1.5 


(u)  =  {[A(e)  A3(ki)~A3]ki}12 

and  the  turbulence  intensity,  Ty,  is  then  calculated  using  the  towing 
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Table  5.1  -  Ocean  Viscous  Dissipation  Rates 

LQ.ca.tiQD  £  r  W/m  Depth  -  m 

Continental  Shelf 

off  Nova  Scotia  10  x  10  to  3  x  10-3  0  to  23. 

Santa  Maria  Azores 

Atlantic  Ocean  3  x  10  6 to  96  x  10” 3  0  to  750. 

Figure  5.1  shows  the  ocean  turbulence  intensity  levels,  Tu,  for  these 
two  locations  and  depth  ranges  as  a  function  of  towing  speed.  The 
resulting  turbulence  intensities,  for  a  towing  speed  of  5  m/s,  vary 
between  0.09  and  0.27  percent  over  a  depth  range  from  0  to  750  m  in 
the  first  location,  while  it  varies  from  0.6  to  8.7  percent  over  a 
depth  range  from  0  to  23  m  in  the  second  location. 

The  wind  tunnel  turbulence  level  could  have  been  reduced  by  the 
addition  of  screens  at  the  inlet;  however,  this  would  have  been  at 
the  expense  of  a  decrease  in  tunnel  maximum  operating  speed  (maximum 
Reynolds  number).  Because  of  the  wide  range  of  variation  in  ocean 
turbulence  intensity,  and  since  maximum  Reynolds  number  in  the  tunnel 
was  considered  more  important,  the  tunnel  was  left  in  its  present 
condition  and  the  wind  tunnel  data  are  presented  as  measured.  If 
desired,  a  correction  for  turbulence  intensity  on  the  skin  friction 
coefficient  is  available  from  Raghunathan  [19] : 

ACf/Cf  =  12.6(T  J1*4  (5.4) 

o  u 

for  0  <  Tu  <  0.0  4 

where  Cf  =  skin  friction  coefficient  at  zero  turbulence. 

o 

5.3.  Wind  Tunnel  Models 

All  five  of  the  fairing  models  were  of  the  same  size  and 
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construction.  Each  model  had  a  chord  length  of  30.5  cm,  a  span  of  [ 

i 

78.11  cm  and  thicknesses  of  7.6  cm  for  the  Fathom  fairing  and  12.2  cm  ; 
for  the  other  four  sections.  The  models  were  made  using  a  urethane  I 

i 

l 

foam  core  with  a  birch  plywood  and  fiberglass  skin.  The  urethane  foam  j 

core  was  cut  to  the  fairing  shape  on  a  numerically  controlled  milling  j 

I 

machine.  The  foam  core  was  then  cut  to  accept  mahogany  supporting  I 

members  that  formed  the  leading  edge,  trailing  edge,  and  side  supports  i 

at  mid-chord.  Profiled  birch  plywood  with  a  recessed  brass  plate 

i 

acted  as  a  mount  for  the  pivot  pin  bearings,  Fig.  5.2.  | 

t 

The  five  models  are  shown  in  Fig.  5.3.  Sixteen  pivot  points 

along  the  chord  at  five  percent  chord  intervals  were  available  for  ! 

1 

each  model.  The  Liebeck  model  was  also  tested  in  two  additional  I 
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configurations,  the  trailing  edge  wedge  and  the  vortex  generator 
strip.  The  trailing  edge  wedge,  Fig.  5.4,  had  a  length  of  0.05  chord 


with  an  included  angle  of  90  degrees  and  a  (t/c)  =  0.105.  The 

Ddi  SG 


vortex  generator  configuration.  Fig.  5.5,  was  based  on  the 
recommendations  of  Pearcey  [20],  where  the  length  1,  height,  h,  and 
spacing,  d,  were  0.01,  0.025  and  0.1  chord.  The  vortex  generators 
were  fabricated  in  the  form  of  a  continuous  strip  which  could  be 
positioned  at  any  station  along  the  chord.  The  tests  were  conducted 
with  the  strip  at  x/c  =  0.35. 

The  two-dimensional  models  were  mounted  in  a  vertical  position 
between  the  top  and  bottom  of  the  wind  tunnel,  Fig.  5.6,  and  were  free 
to  rotate  about  a  pivot  aligned  with  the  leading  edge.  The  pivot 
position  was  variable  along  the  section  chord,  with  the  center  of 
rotation  positions  at  five  percent  chord  stations  along  the  chord. 
Brass  pivot  pins  were  positioned  in  the  top  and  bottom  plates  of  the 
tunnel  test  section  and  were  fitted  into  the  brass  plates  that  were 
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recessed  into  the  top  and  bottom  faces  of  the  models.  A  clear 
plexiglass  plate  acted  as  the  top  viewing  window  for  the  wind  tunnel 
test  section.  Several  holes  and  a  curved  slot  were  cut  in  the 
plexiglass  window  to  hold  the  pivot  pin  and  locking  pin  and  to  mount 
the  moment  arm  for  torque  measurements.  A  protractor  was  etched  into 
the  plexiglass  so  that  the  yaw  angle  of  the  model  could  be  measured. 
The  locking  pin  was  used  to  prevent  the  rotation  of  the  model  and  to 
align  it  with  the  flow  when  the  wake  survey  was  being  conducted. 

5.4.  Data  Acquisition  System 

An  Apple  II  computer  was  used  as  the  data  acquisition  system. 
Fig.  5.7.  This  system  was  chosen  in  part  because  of  the  availability 
of  expansion  hardware  and  software,  and  partly  for  the  ease  with  which 
it  could  be  transported  to  the  non-instrumented  wind  tunnel. 
Table  5.2  lists  the  expansion  cards  that  were  added  to  the  Apple  II  to 
render  it  capable  of  handling  the  acquisition,  reduction,  storage,  and 
output  of  the  data.  The  slot  number  indicates  the  card's  location  on 
the  Apple  bus. 


Table  5.2  -  Apple  II  Computer  Expansion  Hardware 


Apple  XI  Bus  Card  Name 


Manufacturer  Function 


Slot  0 
Slot  1 
Slot  2 
Slot  3 
Slot  4 
Slot  5 
Slot  6 
Slot  7 


Ramcard 
Grappler  + 
SuperSerial  Card 
Videoterm 
AD  213 

Digisector  DS-65 
Disk  II  Interface 
Sof tcard 


Microsoft 

OrangeMicro 

Apple 

Videx 

Tec-Mar 

MicroWorks 

Apple 

Microsoft 


64K  Memory 
Printer  Interface 
Serial  Communictions 
80-Column  Display 
A/D  Motherboard 
Digitizer 

DOS  of  2  Disk  Drives 
Z80A  CPU 


Peripherals  included  a  Zenith  monochrome  video  monitor,  a  Tec-Mar 
12-bit  A/D  daughter  board,  a  C.Itoh  Model  8510A  dot  matrix  printer, 
and  an  HP  7470A  graphics  plotter.  The  A/D  converter  was  used  for  the 


Fig.  5. 
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7.  Apple  II  microcomputer  data  acquisition  system 
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collection  of  the  pressure  and  force  data  during  the  wind  tunnel 
tests.  Pressure  measurements  were  made  for  the  determination  of  the 
fairing  drag  and  for  recording  upstream  dynamic  pressure.  Force  data 
were  taken  to  determine  che  restoring  moment  coefficients.  These 
measurements  were  not  made  simultaneously.  The  wake  pressure  survey 
was  done  with  the  fairing  model  pinned  to  prevent  rotation.  The 
models  were  free  to  pivot  for  the  measurement  of  the  restoring  or 
destabilizing  moments. 

The  analog  output  of  the  pressure  or  force  transducers  was  input 
to  the  A/D  conversion  circuitry.  The  A/D  board  was  manufactured  by 
Tec-Mar,  Inc.  The  signal  conversion  was  handled  with  12-bit  accuracy 
at  a  rate  of  30  kHz.  The  A/D  daughter  board  included  an  Intel  8253 
programmable  interval  timer  (PIT)  for  timing  and  interrupt  control. 
The  analog  input  range  was  +10  volts.  The  digital  output  was  in  twos 
complement  binary  format.  A  Tec-Mar  manufactured  mother  board  handled 
the  interfacing  between  the  A/D  and  the  Apple  bus.  Fig.  5.7  shows  a 
schematic  of  the  A/D  conversion  system. 

Software  for  data  acquisition  was  written  in  FORTRAN.  The 
programs  were  compiled  with  a  Microsoft  compiler  written  for  the 
Apple  II.  Monies  of  the  data  acquisition  and  reduction  programs  are 
contained  in  Appendices  A  and  B.  For  both  pressure  and  force  signals, 
only  the  average  value  at  each  station  was  required.  Thus  the  data 
reduction  programs  did  not  save  each  conversion,  but  summed  the  data 
and  divided  by  the  total  number  of  conversions  to  determine  the 
average  signal.  Five  hundred  conversions  were  made  at  the  maximum 
frequency.  Allowing  for  software  handling,  the  signal  was  averaged 
over  a  period  of  5  seconds.  No  analysis  was  done  on  the  signal 
fluctuations  nor  frequency  components. 
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The  average  digital  signal  was  converted  to  pressure  or  force 
units  using  calibration  constants  that  were  calculated  by  a  separate 
A/C  conversion  program.  The  balance  of  the  data  reduction  programs 
handled  the  operator  interaction  and  interrogation,  the  data 
integration,  the  nondimensionalization,  the  plotting,  and  the  data 
output  functions. 
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6.0.  EXPERIMENTAL  MEASUREMENTS 

The  following  measurements  were  made: 

(1)  boundary  layer  visualization  to  determine  positions  of  laminar  to 
turbulent  transition,  and  laminar  or  turbulent  boundary  layer 
separation 

(2)  wake  survey  to  determine  two-dimensional  section  drag  coefficient 

(3)  determination  of  position  of  hydrodynamic  center,  or  neutral 
stability  point 

(4)  torque  about  the  center  of  rotation  as  a  function  of  yaw  angle. 

6.1.  Boundary  Layer  Visualization 

Two  techniques  were  used  for  flow  visualization,  including 
fluorescent  mini-tufts  and  an  oil/ultraviolet  light  method.  The 
fluorescent  mini-tuft  technique,  developed  by  Crowder  [21],  utilizes 
extremely  thin  fluorescent  nylon  monofilaments,  about  0.018  mm  in 
diameter,  attached  to  the  model  with  a  drop  of  lacquer  adhesive.  The 
use  of  the  mini-tufts  results  in  negligible  interference  with  the 
flow,  as  determined  by  wind  tunnel  tests  [21,22].  Thus  a  large  matrix 
of  tufts  may  be  used  for  flow  visualization  without  changing  the  flow 
patterns.  The  tufts  were  used  to  map  ♦•’-e  areas  of  turbulent  boundary 
layer  separation.  The  tufts,  approxiiu-cely  2.5  cm  in  length,  must  be 
viewed  under  ultraviolet  light  to  contrast  them  with  the  white  model 
surface.  A  matrix  of  the  tufts  was  applied  at  10  percent  chord 
stations  along  the  chord,  with  approximately  equal  spacings  along  the 
span  of  the  model.  With  a  white  model  background  viewed  under 
ultraviolet  light,  the  model  surface  disappeared  and  the  tufts  were 
illuminated. 


The  second  technique,  which  involved  the  use  of  oil  and 
ultraviolet  light,  was  used  to  determine  the  position  of  transition 


from  laminar  to  turbulent  flow  in  the  boundary  layer,  the  existence  of 
a  laminar  separation  bubble  and  turbulent  separation.  The  latter  thus 
serves  as  a  check  on  the  fluorescent  mini-tuft  method.  A  mixture  of 
kerosene  (70  percent  by  volume),  104W0  motor  oil  (20  percent)  and 
vinegar  (10  percent)  was  used  along  with  approximately  four  to  five 
tablespoons  of  fluorescent  yellow  pigment  per  pint  of  the  liquid.  The  ( 

kerosene  provides  the  base  for  the  mixture,  the  motor  oil  slows  down 

the  rate  of  evaporation  of  the  kerosene  in  the  wind  stream,  and  the 
vinegar  is  used  to  keep  the  powdered  fluorescent  pigment  in  | 

suspension.  It  has  been  found  through  experiments  that  a  lemon  yellow 
pigment  used  in  combination  with  a  white  model  background  results  in 
the  highest  contrast.  I 

Both  types  of  flow  visualization  were  viewed  under  100-watt 
ultraviolet  lights  approximately  0.5  to  1.0  m  from  the  model. 
Photography  was  accomplished  with  settings  of  f/2  to  f/4  at  1/15  to 
1/30  of  a  second,  using  a  Minolta  ultraviolet  light  filter  over  a 
Griffin  1-A  skylight  filter. 

6.2.  Drag  Coefficient 

Two  wind  tunnel  methodologies  are  available  for  determining  the 
drag  coefficient  of  streamline  shapes.  One  technique  requires  the  use 
of  mechanical  or  electrical  balances  to  which  the  model  is  attached, 
for  direct  measurement  of  the  lift,  drag,  and  pitching  moment 
aerodynamic  forces.  The  other  technique  is  indirect.  The  profile 
drag  is  determined  from  the  velocity  distribution  in  the  wake  of  the 
fairing.  This  method  was  chosen  for  this  study  because  it  was 
believed  that  the  position  of  the  hydrodynamic  center  could  be 
determined  with  greater  accuracy  using  the  free  pivoting  technique. 

This  precluded  the  use  of  models  that  would  fit  on  the  available  wind 
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tunnel  balance. 


The  momentum  loss  method  for  determining  the  profile  drag  is  well 
known,  Pope  [15].  When  the  momentum  equation  is  integrated  around  a 
control  volume  that  includes  the  wake,  it  can  be  shown  that  the  drag 
force  on  the  body  is  given  by  the  expression: 

oo 

D  =  bp  /  U(Uoo  -  U)dy  .  (6.1) 

—  OO 

6.2.1.  Momentum  Wake  Rake 

This  integration  must  be  done  with  the  wake  profile  data  taken 
far  enough  downstream  from  the  body  so  that  the  static  pressure  at  the 
measuring  section  is  equal  to  that  in  the  undisturbed  stream. 

Pope  [15]  states  that  locating  the  wake  rake  at  least  0.7  chord  behind 

the  trailing  edge  of  the  wing  is  sufficient. 

The  wake  rake,  shown  in  Fig.  6.1,  was  positioned  at  mid-tunnel 
height  and  0.678  chord,  20.3  cm,  length  downstream  of  the  model 
trailing  edge.  The  wake  rake  was  fixed  in  this  position  and  was  not 

adjustable  to  port  or  starboard.  Thirty-one  total  head  probes  were 

spaced  at  0.5  in.  (1.27  cm)  intervals  along  the  wake  rake.  The 

center  5  probes  were  spaced  38.1  cm  apart.  The  total  width  of  the 

wake  survey  was  38.1  cm.  The  probes,  0.16  cm  o.d.  rigid  stainless 
steel  tubing,  extended  7.62  cm  from  the  body  of  the  rake.  Static 

pressure  probes  were  made  from  the  same  tubing.  These  tubes  extended 

10.16  cm  from  the  rake  body  and  the  ends  were  sealed  and  rounded. 
Four  0.08  cm  holes  were  drilled  perpendicular  to  the  axis  of  the 
tubing,  2.54  cm  downstream  from  the  leading  edge.  The  static  probes 
were  spaced  at  5.08  cm  intervals  over  the  center  31.75  cm  of  the  wake 
rake.  The  static  pressures  were  taken  in  the  same  plane  and  the  same 
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distance  downstream  from  the  model  as  the  total  head  measurements. 

Flexible,  multi-colored  "strip-a-tube"  was  lead  through  the  body 
of  the  wake  rake,  down  through  the  faired  supporting  struts,  and  out 
to  a  distribution  panel.  Each  pressure  tap  of  the  wake  rake  was 
individually  connected  by  way  of  a  valve  system  to  a  differential 
pressure  transducer.  Six  of  the  ten  static  pressure  probes  from  the 
central  region  of  the  wake  were  attacked  to  one  side  of  the 
transducer,  thus  providing  the  average  static  pressure  in  the  wake. 
The  total  head  probes  were  then  individually  input  to  the  other  side 
of  the  transducer.  The  output  signal  was  the  dynanic  pressure  at  each 
station  in  the  wake,  (q  =  H  -  Pg) . 

The  transducer  was  a  Dynasciences  model  P109D  0  to  ±15  PSID 

variable  reluctance  type  differential  pressure  transducer.  The 
transducer  was  driven  by  a  Dynasciences  model  CD10  DC  output 
Carrier-Demodulator  operating  at  10  DC  volts  full  scale.  The  wake 
dynamic  pressures  varied  from  0.003  PSIG  to  0.1  PSIG.  The  output 
voltage  ranged  from  1.0  to  3.0  volts.  Calibration  curves  were 
calculated  before  each  set  of  wake  surveys  using  a  kerosene  manometer. 
The  calibrations  were  linear  to  ±0.1  percent  of  best  straight  line 
fit.  The  output  voltages  from  the  transducer  were  input  to  the  A/D 
circuit  described  in  Section  5.4. 

The  data  reduction  program  prompted  the  operator  to  open  each 
probe  valve  in  sequence,  and  then  collected,  converted,  and  stored  the 
pressure  data.  The  data  collected  from  the  port  and  starboard  side  of 
the  wake  were  folded  and  averaged  so  that  the  wake  was  forced  to  be 
symmetrical.  The  integration  required  by  Eqn.  6.9  was  accomplished 
using  the  Fr itsch-Carlson  [23]  formulas  for  piecewise  cubic  data 


interpolation 


and  Gauss-Legendr e  Quadrature 


numerical 


,  i*.  j 


integration  [ 24 ] .  The  limits  of  the  integration  (the  wake  width)  were 
determined  by  operator  inspection  of  the  wake  profile  shape.  The 
program  plotted  the  wake  profile  and  the  operator  would  then  choose 
and  input  the  integration  limits.  A  sample  wake  survey  profile  for 

the  Liebeck  section  is  shown  in  Fig.  6.2  for  r  =  0.511  x  106 .  The 

c 

symmetry  of  the  wake  profile  was  due  to  the  folding  and  averaging 
described  above. 

The  experimental  parameters  (model  name  and  configuration,  date, 
model  dimensions,  rake  position,  and  probe  locations)  were  read  into 
the  program  from  a  special  parameter  file.  The  wind  tunnel  operating 
conditions  of  temperature  and  upstream  dynamic  pressure  were  input  by 
the  program  operator.  The  program  then  calculated  the  downstream  R_, 

C 

the  blockage  correction  factor  used  in  Eqn.  6.9,  and  the  drag 
coefficient. 

The  drag  coefficient  may  be  based  on  planform  area  (be)  as  is 
common  in  aeronautical  work,  or  on  frontal  area  (bt),  as  is  common 
practice  among  marine  engineers.  Thus 


=  DApirbc 


(6.2) 


=  DA  U2bt 


(6.3) 


Throughout  this  paper  the  drag  coefficients  will  be  given  based  on 
frontal  area.  Using  Equation  6.2,  the  drag  coefficient  equation 
becomes: 


'd  ..  . 


~  t:  /  C  (q/qj  1/2  -  (q/qj]dy 


(6.4) 


UlflKE  POSITION  (  in.  ) 


Fig.  6. 


2.  Wake  profile  for  Liebeck  section  at  R  =  0.511  x  10 


6.2.2.  Blockage  Correction 


The  flow  conditions  that  the  model  experiences  in  a  wind  tunnel 
are  not  the  same  as  those  in  free  air  or  the  open  ocean.  The  presence 
of  the  wind  tunnel  walls  changes  the  velocity  of  air  flow  over  the 
model.  The  two  corrections  that  are  important  in  the  calculation  of 
drag  coefficients  are  "solid  blockage"  and  "wake  blockage."  Pope  [15] 
defines  the  total  blockage  as 


E 


E  .  +  E  , 

sb  wb 


(6.5) 


where  E  is  the  sum  of  the  solid  blockage,  ESb'  an<*  wa*te  blockage, 
Ewb-  Or,  in  terms  of  the  change  in  freestream  velocity  due  to  the 
presence  of  the  model 


E  =  (U0/Ucc)  ■  1 


(6.6) 


or 

E  =  ^%7qT-  1  (6.7) 

where  UQ  and  qQ  are  the  velocity  and  dynamic  pressure  downstream  of 
the  model  at  the  edge  of  the  wake. 

The  correction  of  the  wake  survey  drag  coefficient  for  blockage 
by  the  method  proposed  by  Shaw,  Sotos  and  Solano  [23]  of  NACA  Lewis 
Research  Center.  They  employed  the  following  procedure  to  correct  for 
the  effects  of  blockage. 

(1)  The  original  wake  survey  is  examined  to  determine  the  endpoints 
of  the  viscous  wake. 

(2)  The  values  of  the  velocity  ratio  (l^/U^)  for  the  two  wake 
endpoints  are  compared,  and  the  larger  volume  used  to  calculate  a 
velocity  ratio  correction  factor  by  subtracting  1.0  from  this 


value. 


I 


(3)  All  velocity  ratio  values  contained  within  the  viscous  wake  are 
adjusted  by  subtracting  the  velocity  ratio  correction  factor 
prior  to  computing  the  section  drag  coefficient. 

The  wake  data  were  corrected  for  blockage  in  the  following 
manner.  The  blockage  correction  factor  was  first  determined 
experimentally  as 

E  =  (Uo/Uj  -  1  =  AU/U^  (6.8) 


where:  =  upstream  velocity 

Uq  =  downstream  velocity  at  edge  of  wake. 

The  drag  coefficient  was  then  computed  using  the  NASA  method 
described  in  Shaw,  et  al.  [24],  which  includes  the  blockage  correction 
factor . 


'  \ 

d(bo) 

m 

where:  U 

* 

dy 

=  §  /  [(U/U^-E)  (1  -~  +  E)]dy 


dy  =  distance  across  wake. 


(6.9) 


6.2.3.  Base  Drag  Coefficient 

Static  pressures  at  the  blunt  trailing  edge  were  measured  for 
three  of  the  fairing  section  shapes.  A  static  pressure  probe  was 
mounted  along  the  centerline  of  the  trailing  edge.  Pressure  taps  were 
drilled  in  this  probe  at  six-inch  intervals  over  the  center  18  inches 
of  the  fairing  span.  The  static  pressures  were  measured  relative  to 
the  free  stream  static  pressures  which  were  measured  downstream  from 
the  model  at  the  wake  rake. 

The  contribution  to  the  total  drag  of  the  blunt  trailing  edge  was 
calculated  from  the  following  equations. 


PsA/qbc 


(6.10) 


where: 

C.  =  base  drag  coefficient, 

ate 

P_  =  base  static  pressure 

A  =  bt  =  base  area  of  trailing  edge, 
q  =  downstream  dynamic  pressure  corrected  for 

blockage  effects, 
be  =  fairing  planform  area. 

The  ratio  C,  /C.  is  the  fraction  of  the  total  drag  due  to  the 

ate  (be) 

blunt  trailing  edge. 

6.3.  Hydrodynamic  Center 
6.3.1.  Definition 

The  aerodynamic  center  is  defined  as  the  position  along  the 
airfoil  chord  for  which  the  moment  coefficient  ,  Cm=„,  is  constant. 


that  is,  independent  of  CT. 


a.c. 


=  n  -  dC  /dCT 
m  L 
n 


(6.11) 


where 


Cm  =  constant  (  =  0.  for  symmetrical  section) 

a.c. 

n  =  moment  center  location  along  chord 

O  =  moment  coefficient  about  n 

m 

n 

C  =  M  /qSC 

n 

CL  *  lift  coefficient. 

For  symmetrical  sections,  the  value  of  the  pitching  moment 

coefficient,  C  ,  is  zero. 

m 

n 

The  neutral  point  is  defined  as  that  point  along  the  chord  where 
the  slope  of  the  pitching  moment  versus  lift  coefficient  curve  is 
equal  to  zero,  i.e., 


n  =  neutral  point  location. 

Thus  for  a  symmetrical  section, 

a.c.  =  n  -  dC  /dC  (6.13) 

mn  L 

=  n  -  0 

and 

ci  •  c  •  n  • 

Therefore,  for  symmetrical  sections,  the  neutral  point  and  the 
aerodynamic  center  are  one  and  the  same.  In  this  study  the  terms 
aerodynamic  center  and  hydrodynamic  center  are  used  interchangeably. 

6.3.2.  Free  Pivot  Technique 

Two  methods  were  used  to  find  the  chordwise  position  of  the 

hydrodynamic  center.  In  the  first  method,  a  visual  measurement  was 
made  of  the  equilibrium  angle  of  attack  of  the  fairing  as  a  function 
of  the  mechanical  pivot  location.  When  the  pivot  point  is  forward  of 
the  hydrodynamic  center,  the  fairing  will  remain  aligned  with  the 

flow,  exhibiting  "weathervane"  stability  and  the  angle  of  attack  will 
be  zero  degrees. 

The  pivot  location  was  moved  progressively  aft  in  steps  of 

0.05  c.  At  each  pivot  position  the  equilibrium  angle  of  attack  was 
recorded.  The  angles  were  measured  using  the  protractor  that  was 

inscribed  in  the  plexiglass  top  of  the  wind  tunnel.  The  accuracy  of 
the  readings  was  ±0.5  deg.  When  the  pivot  position  is  moved  aft  of 

the  hydrodynamic  center,  the  fairing  will  assume  an  equilibrium  angle 
of  attack  other  than  zero  degrees.  This  equilibrium  angle  will  be 
positive  or  negative,  depending  on  the  direction  of  the  initial 


disturbance.  Because  the  step  size  between  pivot  positions  was 
0.05  c,  the  hydrodynamic  center  was  located  by  this  technique  with  an 
accuracy  of  +0.025  c. 

6.3.3.  Pivot  Point  Torque  Technique 

Restoring  or  destabilizing  moments  were  measured  about  several 
pivot  points  as  an  alternative  method  of  locating  the  hydrodynamic 
center.  When  the  mechanical  center  is  forward  of  the  hydrodynamic 
center,  the  moments  tend  to  align  to  fairing  with  the  flow.  When  the 
mechanical  center  is  aft  of  the  aerodynamic  center,  the  moments  are 
destabilizing  and  cause  misalignment  with  the  flow.  Thus  the  sign  of 
the  moment  coefficient  changes  from  negative  to  positive  as  the  pivot 
location  moves  aft  past  the  hydrodynamic  center. 

The  moments  were  measured  using  a  free  pivot  technique.  The 
fairing  model,  free  to  to  pivot  about  a  selected  chordwise  position, 
was  rotated  to  an  equilibrium  angle  of  attack  to  the  flow  with  a  force 
transducer  to  measure  the  restoring  or  destabilizing  moment  about  the 
pivot  location. 

The  output  signal  from  the  force  transducer  was  connected  to  an 
A/D  converter  driven  by  the  Apple  computer.  Moment  data  were  taken 
for  each  fairing,  at  pivot  locations  ranging  from  5  to  35  percent  of 
the  chord  length,  and  for  geometric  angles  of  attack  ranging  0  to 
20  degrees  in  increments  of  two  degrees.  The  upstream  dynamic 
pressure  was  approximately  8.0  psf  (Reynolds  number  0.5  x  106). 

The  data  reduction  program  calculates  and  prints  the  moment 
coefficient  as  a  function  of  angle  of  attack  for  each  pivot  location. 
The  moment  coefficient  is 


C 

m 


Fr/qbc2 


(6.14) 


where  r  is  the  moment  arm  and  F  is  the  measured  force. 

Because  of  the  presence  of  the  wind  tunnel  walls,  the  measured 
moment  coefficients  differ  from  the  moments  that  would  be  experienced 
by  a  fairing  in  the  open  ocean.  A  correction  is  necessary  due  to  the 
effects  of  blockage  and  streamline  curvature  that  results  from  the 
presence  of  the  walls.  The  measured  moment  coefficients  are  greater 
than  they  would  be  in  the  freestream  environment. 

Pope  [15]  gives  the  moment  coefficient  correction  as 

C  =  C  (1  -  2E)  +  oC./4  (6.15) 

m  m  1 

n 

with 

o  =  ( 7T 2 / 4 8 )  (C  (h)  2  =  0.04. 

The  subscript,  u,  indicates  the  uncorrected  coefficient.  The 
correction  to  the  moment  coefficient  is  a  function  of  blockage,  E,  and 
the  lift  coefficient,  Cj_.  The  lift  coefficient  is  a  function  of  angle 
of  attack,  a.  Likewise,  the  measured  angle  of  attack  must  be 
corrected  for  the  effects  of  streamline  curvature.  Pope  gives  this 
correction  as 

a  =  a  +  +  4C  ).  (6.16: 

u  2lT  Xu  mn 

With  the  free  pivot  technique  used  in  this  project,  it  was  not 
possible  to  measure  the  lift  coefficients  and  the  blockage  factor  upon 
which  both  these  corrections  depend.  Thus  the  presented  moment  data 
are  uncorrected  for  wind  tunnel  boundary  conditions  and  can  only  be 
used  as  indications  of  the  location  of  the  hydrodynamic  center.  For 
small  angles  of  attack,  less  that  ±60  deg.,  the  moment  data  have  been 


corrected  for  the  effects  of  blockage.  For  these  small  angles,  was 
taken  from  the  drag  measurement  data  and  C  was  assumed  to  be  near 
zero.  Thus  Equation  6.15  becomes 


Cm  =  Cm  -  2V 
u 


(6.17) 


where  EQ  is  the  blockage  measured  at  zero  angle  of  attack. 

Corrections  to  the  measured  angles  of  attack,  au,  were  not 
calculated.  The  resulting  error  in  ranges  from  0  to  5  percent  as 
ranges  from  0  to  20  deg.,  assuming  a  maximum  lift  coefficient  of  1.0. 
and  maximum  moment  coefficient  of  0.5. 


7.0.  EXPERIMENTAL  RESULTS 

7.1.  Liebeck  Section 

7.1.1.  Boundary  Layer 

The  results  of  both  boundary  layer  visualization  surveys  for  all 
three  Liebeck  model  configurations  are  shown  in  Figs.  7.1,  7.2.,  and 
7.3  as  a  function  of  Reynolds  number.  Noted  are  the  chordwise 
position  of  the  transition  from  laminar  to  turbulent  boundary  layer, 
the  existence  of  a  laminar  separation  bubble  and  turbulent  boundary 
layer  separation.  A  laminar  separation  bubble  was  found  to  exist  at 
the  transition  from  laminar  to  turbulent  flow.  The  position  of  this 
laminar  separation  bubble  coincided  approximately  with  the  point  of 
minimum  pressure  coefficient,  which  is  at  approximately  40  percent  of 
the  chord.  It  is  seen  from  Fig.  7.1  that  past  a  Reynolds  number  of 
4.0  x  105,  the  turbulent  separation  over  the  aft  portion  of  the 

section  has  completely  disappeared.  At  the  lower  Reynolds  numbers, 

approximately  1  x  105  to  2.7  x  105,  the  aft  portion  was  separated  from 
about  40  percent  chord  on.  The  boundary  layer  studies  were  conducted 
simultaneously  on  either  side  of  the  model. 

The  vortex  generators,  which  were  fixed  at  0.35  c,  eliminated  the 
turbulent  separation  at  all  but  the  lowest  Reynolds  number,  Fig.  7.2. 
The  addition  of  the  trailing  edge  wedge,  Fig.  7.3,  had  no  effect 
except  at  the  higher  Reynolds  numbers.  A  standing  vortex  flow  was 
generated  in  the  region  of  the  wedge,  such  that  a  turbulent  boundary 
separation  condition  was  initiated.  In  all  three  cases,  the 
transition  from  laminar  to  turbulent  flow  occurred  at  the  same 
location.  Transition  occurred  at  0.30  c  for  the  low  Reynolds  number 
of  1.2  xlO5  and  moved  progressively  aft  to  0.47  c  at  5.5  x  105.  The 

effect  of  the  massive  turbulent  separation  from  1.0  to 
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2.7  x  105  apparently  caused  this  migration  forward. 

The  results  of  the  fluorescent  mini-tuft  survey.  Figs.  7.4  and 
7.5,  show  the  fully  attached  flow  for  the  baseline  configuration  at  a 
Reynolds  number  of  3.9  x  10 5 ,  and  the  effect  of  the  standing  vortex  on 
the  trailing  edge  wedge  flow  pattern.  The  separated  flow  condition 
which  existed  for  a  Reynolds  number  of  1.5  x  105  is  shown  in  Fig.  7.6. 
Separation  starts  at  approximately  the  45  percent  chord  position, 
after  which  the  motion  of  the  mini-tufts  is  quite  dramatic.  The  oil 
technique  results  are  shown  in  Fig.  7.7  for  a  Reynolds  number  of 
4.0  x  105.  The  beginning  of  the  laminar  separation  bubble  can  be  seen 
at  the  40  percent  chord  station.  The  fully  attached  turbulent 
boundary  layer  aft  of  this  location  can  be  distinguished  from  the 
laminar  boundary  layer  forward  of  this  location  by  the  relative 
contrast  between  the  pigment  coating  which  remains  on  the  surface. 
Since  the  shearing  stress  levels  are  lower  in  the  laminar  boundary 
layer,  less  of  the  oil  is  scrubbed  away,  and  consequently  the  yellow 
pigment  has  a  much  higher  contrast  in  this  area.  In  the  turbulent 
boundary  layer,  the  shearing  stress  levels  are  higher  and  consequently 
more  of  the  liquid  is  scrubbed  away,  resulting  in  a  darker  area. 


7.1.2.  Drag  Coefficient 

The  blockage  correction  factors  are  shown  in  Fig.  7.8  as  a 
function  of  Reynolds  number  based  on  UQ.  Blockage  over  the  Reynolds 
number  range  tested  varies  from  57  down  to  2  percent  for  the  baseline 
model,  while  it  was  constant  at  about  13  percent  for  the  trailing  edge 
wedge  configuration.  The  addition  of  the  vortex  generators  to  the 
baseline  model  decreased  the  blockage  due  to  modification  of  the 
boundary  layer.  A  typical  baseline  model  wake  profile  for  a  Reynolds 
number  of  0.511  x  105  is  shown  in  Fig.  6.2.  The  symmetry  of  the 
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profile  is  due  to  port  and  starboard  averaging  of  the  data. 

The  two-dimensional  drag  coefficient  corrected  for  blockage  for 
the  three  configurations  tested  is  shown  in  Fig.  7.9  as  a  function  of 
the  Reynolds  number  based  on  chord.  In  order  to  provide  a  comparison, 
the  results  from  Hoerner  [11]  are  included  for  symmetrical  sections 
with  thickness/chord  ratios  of  37  and  50  percent.  It  should  be  noted 
that  the  drag  coefficient  is  based  on  the  projected  frontal  area 
rather  than  the  planform  area.  The  shape  of  the  baseline  section 
curve  is  seen  to  agree  well  with  the  Hoerner  curve.  The  drag 
coefficient  decreases  down  to  a  Reynolds  number  of  4  x  105,  after 
which  it  is  approximately  constant. 

The  addition  of  the  vortex  generators,  which  eliminated  the 
boundary  layer  separation  at  the  low  Reynolds  number,  is  shown  to 
decrease  the  drag  coefficient  in  this  region.  The  crossover  point 
occurs  at  a  Reynolds  number  of  about  3.0  x  10 5 ,  which  agrees  with  the 
results  shown  in  Fig.  7.1,  where  the  turbulent  separation  disappeared 
on  the  baseline  model.  At  Reynolds  numbers  above  this  point,  the 
model  with  the  vortex  generators  had  a  higher  drag  coefficient, 
indicating  that  the  drag  of  the  vortex  generators  themselves  was  quite 
large  when  working  in  a  range  where  they  had  no  effect  on  the  boundary 
layer. 

The  addition  of  the  trailing  edge  wedge  is  seen  to  increase  the 
drag  coefficient  over  most  of  the  Reynolds  number  range.  At  the 
higher  Reynolds  numbers,  the  drag  coefficient  is  higher  by  a  factor  of 
three. 

For  the  Liebeck  40  percent  section  with  the  trailing  edge  wedge, 
the  percentage  of  trailing  edge  drag,  Fig.  7.10,  varied  from  15  to 
56  percent  for  the  Reynolds  numbers  tested.  For  this  case,  the 
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contribution  to  the  total  drag  of  the  trailing  edge  wedge  was  not 
constant  but  varied  significantly  with  Reynolds  number.  At  the  low 
Reynolds  number  of  0.13  x  10%  the  static  pressure  was  greater  than 
the  freestream  static  pressure,  indicating  separated  flow  at  the 
trailing  edge.  The  boundary  layer  visualization  indicated  that 
separation  had  occurred  for  all  Reynolds  numbers  less  than  0.35  x  106. 

7.1.3.  Hydrodynamic  Center 

The  results  of  the  hydrodynamic  center  (neutral  point)  location 
study  are  shown  in  Fig.  7.11  for  a  Reynolds  number  of  4.8  x  105.  As 
can  be  seen  for  the  baseline  section,  the  hydrodynamic  center  is  quite 
far  forward,  at  about  20  percent  of  the  chord.  For  pivot  positions 
aft  of  this  station,  the  model  would  assume  a  bistatically  stable 
position  at  either  a  plus  or  a  minus  yaw  angle,  depending  on  which  way 
the  initial  disturbance  was  directed.  The  unstable  region  increases 
until  finally,  past  stations  at  approximately  30  percent  of  the 
section  chord,  it  widely  diverges. 

The  effect  of  adding  the  trailing  edge  wedge,  Fig.  7.11,  is 
dramatic.  The  model  is  seen  to  be  stable  for  pivot  positions  all  the 
way  up  to  35  percent  of  the  chord,  past  which  point  it  was 
bistatically  unstable.  This  indicates  that  the  addition  of  the 
trailing  edge  wedge  would  allow  the  use  of  a  larger  diameter  cylinder. 
The  penalty  to  be  paid  for  this,  of  course,  is  the  higher  drag 
coefficient  over  the  operating  Reynolds  number  range. 

The  restoring  torque  coefficient  is  shown  in  Fig.  7.12  as  a 
function  of  pivot  location  along  the  chord  length  for  each  yaw  angle. 
The  point  where  each  curve  crosses  the  zero  moment  axis  corresponds  to 
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7.11.  Yaw  equilibrium  positions  for  the  Liebeck  sections 
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7.2.  Fathom  Fairing  Section 
7.2.1.  Boundary  Layer 


The  boundary  layer  survey  results  for  the  Fathom  section  shown  in 
Fig.  7.13.  Transition  occurred  at  the  20  percent  chord  position  and 
moved  forward  with  increasing  Reynolds  number  to  the  15  percent  chord 
position.  For  this  section  the  laminar  to  turbulent  boundary  layer 
transition  occurred  downstream  of  the  minimum  pressure  point,  which 
was  located  at  the  9. 8  percent  chord  position.  A  laminar  separation 
bubble  also  occurred  at  the  transition  location.  The  bubble  length, 
approximately  10  percent  chord  length  at  the  lowest  Reynolds  number, 
decreased  in  length  as  Reynolds  number  increased. 

Turbulent  separation  occurred  just  before  the  trailing  edge  at 
the  95  percent  chord  position  for  all  Reynolds  numbers  tested.  This 
was  in  agreement  with  the  results  of  Henderson  [2],  who  also  found 
separation  at  the  95  percent  chord  position  at  a  Reynolds  number  of 
0.22  x  106 . 

7.2.2.  Drag  Coefficient 

The  measured  blockage  correction  factors  for  the  Fathom  Flexnose 
section  are  shown  in  Fig.  7.14.  Blockage  ranges  from  a  high  of 
40  percent  at  the  lowest  Reynolds  number  to  12  percent  at  the  highest. 
The  curve  is  quantitatively  similar  to  the  blockage  found  for  the 
Liebeck  sections.  However,  the  blockage  effect  seems  quite  large 
considering  the  section  t/c  ratio. 

Because  there  were  no  major  changes  in  the  behavior  of  the 
boundary  layer  throughout  the  Reynolds  number  range  tested,  the  drag 
coeffients.  Fig.  7.15,  remained  fairly  constant.  The  measured  drag 
coefficients  are  in  good  agreement  with  the  one  data  point  measured  by 
Henderson. 
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Fig.  7.13.  Results  of  the  boundary  layer  survey  of  the  Fathom 
Flexnose  section. 
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7.14.  Measured  blockage  of  the  Fathom  Flexnose  sect' 
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7.2.3.  Hydrodynamic  Center 

Figures  7.16  and  7.17  show  the  results  of  the  stability 
experiment  at  two  different  Reynolds  numbers  for  the  Fathom  fairing, 
0.38  and  0.53  x  106.  The  hydrodynamic  center  is  seen  to  be  located  at 
15  percent  chord  at  the  low  Reynolds  number  and  70  percent  chord  for 
the  high  one.  Regions  of  dynamic  instability  occurred  at  certain 
pivot  point  locations  at  the  low  Reynolds  numbers. 

Moment  coefficient  measurements  plotted  in  Fig.  7.18  were  made  at 
a  Reynolds  nubmer  of  0.495  x  106.  The  hydrodynamic  center  was  found 
to  be  located  near  the  21.5  percent  chord  position,  which  agrees  with 
the  free  pivot  results.  The  hydrodynamic  center  for  the  Fathom 
Flexnose  fairing  was  found  by  Henderson  from  lift  and  moment  data  to 
be  located  at  the  15.3  percent  chord  position  at  a  Reynolds  number  of 
0.22  x  106. 

The  design  mechanical  pivot  position  for  the  prototype  Fathom 
Flexnose  25  percent  fairing  is  the  12  percent  chord  position.  The 
fact  that  the  hydrodynamic  center  is  located  so  closely  behind  the 
mechanical  center  can  explain  most  of  the  unacceptable  behavior  of 
this  fairing  in  sea  trials.  A  complete  description  of  the  problems 
encountered  with  this  fairing  is  given  by  Henderson  [2], 

7.3.  NACA  0040  Section 
7.3.1.  Boundary  Layer 

The  results  of  the  boundary  layer  survey  of  the  NACA  0040  fairing 
are  shown  in  Fig.  7.19.  The  laminar  to  turbulent  boundary  layer 
transition  location  was  found  to  move  forward  toward  the  leading  edge 
with  increasing  Reynolds  numbers  between  0.2  and  0.3  x  106.  At  a 
Reynolds  number  of  0.3  x  106,  transition  occurs  at  the  30  percent 
chord  position.  The  minimum  pressure  point  determined  from  the 
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potential  flow  calculations  is  the  19  percent  chord  location.  The 
maximum  thickness  for  the  NACA  0040  section  is  at  the  30  percent  chord 
location.  Thus  transition  occurred  aft  of  the  minimum  pressure  point 
and  just  at  the  maximum  thickness  position.  For  Reynolds  numbers 
greater  than  0.3  x  106 ,  the  location  of  the  boundary  layer  transition 
was  not  clearly  indicated  with  the  oil/dye  technique.  There  was  no 
laminar  separation  bubble  visible  and  no  apparent  transition  point. 
It  was  possible  that  transition  occurred  at  the  leading  edge.  This 
investigation  was  unable  to  determine  its  location  with  any  certainty. 

Transition  locations  as  a  function  of  Reynolds  number  are  plotted 
in  Fig.  7.20  for  other  NACA  00XX  sections  [8].  In  each  case,  the 
boundary  layer  transition  location  moved  forward  toward  the  leading 
edge  with  increasing  Reynolds  number.  There  was  no  indication  that 
transition  might  occur  at  the  leading  edge  on  any  of  the  thinner  NACA 
sections. 

The  turbulent  boundary  layer  separation  location  varied  slightly 
with  Reynolds  number,  moving  toward  the  trailing  edge  as  the  Reynolds 
numbers  increased.  Throughout  the  Reynolds  number  range  tested,  from 
0.20  to  0.61  x  10 6,  the  boundary  layer  was  never  fully  attached. 

7.3.2.  Drag  Coefficient 

The  measured  blockage  correction  factors  for  the  NACA  0040 
section  are  shown  in  Fig.  7.21.  The  blockage  ranges  from  26.5  percent 
at  the  low  Reynolds  number  of  0.15  x  10 6  to  6.5  percent  at  Reynolds 
number  0.56  x  10 6.  The  results  of  the  two-dimensional  drag 
coefficient  measurements  for  the  NACA  0040  fairing  are  shown  in 
Fig  7.22.  Also  shown  are  drag  coefficients  for  a  Joukowski  40  percent 
section,  and  a  NACA  0033  measured  by  Fage,  Falkner,  and  Walker  [24] 
and  Althaus  [13] . 
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The  drag  of  the  NACA  0040  fairing  decreased  slightly  over  the 
Reynolds  number  range  considered.  The  coefficient  had  a  maximum  value 
of  0.108  at  a  Reynolds  number  of  0.14  x  106  and  a  low  of  0.08  at  a 
Reynolds  number  of  0.56  x  106.  The  slight  decrease  in  drag 

coefficient  with  increasing  Reynolds  number  was  related  to  both  the 
decreasing  portion  of  laminar  boundary  layer  drag  increase,  and 
decreasing  portion  of  separated  turbulent  boundary  layer  drag 
decrease.  These  two  changes  have  opposite  effects  on  the  drag 

coefficient. 

7.3.3.  Hydrodynamic  Center 

Fig.  7.23  shows  the  results  of  the  stability  measurements  for 
this  section  at  a  Reynolds  number  of  0.48  x  106 .  The  location  of  the 
hydrodynamic  center  was  found  to  be  forward  of  the  leading  edge.  With 
the  pivot  location  at  5  percent  of  the  chord  length,  the  yaw 
equilibrium  position  was  +6  deg.  The  yaw  equilibrium  angle  increased 
as  the  pivot  location  was  moved  aft.  With  the  pivot  location  at 
35  percent  chord  the  fairing  was  in  equilibrium  at  a  yaw  angle  of 
+41  deg. 

Measurements  of  the  moment  coefficient,  Fig.  7.24,  also  confirmed 
that  the  location  of  the  hydrodynamic  center  was  forward  of  the 
leading  edge.  Measured  moments  were  positive  for  small  angles  of 
attack,  up  to  10  deg.,  for  pivot  positions  at  Xp/c  =  0.05,  0.10,  and 

0.15.  Positive  moments  are  destabilizing,  resulting  in  misalignment 
of  the  fairing  with  the  flow. 

Fig.  7.25  shows  a  plot  of  the  location  of  the  aerodynamic  center 
vs.  t/c  ratio  for  NACA  00XX  sections  as  given  by  Eastman  [9].  The 
movement  toward  the  leading  edge  with  increasing  t/c  ratio  is  clearly 
indicated.  For  the  NACA  0035  section  the  aerodynamic  center  was  found 
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by  Bullivant  [8]  to  be  at  a.c.  =  0.16.  The  location  of  the 
hydrodynamic  center  forward  of  the  leading  edge  for  the  NACA  0040  was 
a  significant  phenomenon  which  could  not  have  been  predicted  from  the 
investigation  of  the  other  NACA  00XX  sections. 

The  misalignment  or  lack  of  weathervane  stability  of  the 
NACA  0040  section  at  all  mechanical  pivot  locations  precludes  the 
possibility  of  using  this  section  as  a  cable  fairing.  The  large 
hydrodynamic  lift  and  increase  in  drag  due  to  the  equilibrium  angle  of 
attack  would  cause  significant  problems.  The  addition  of  a  trailing 
edge  wedge  as  was  fitted  on  the  Liebeck  40  percent  fairing  would  move 
the  location  of  the  hydrodynamic  center  aft  to  a  more  favorable 
location.  No  measurements  were  made  to  find  the  hydrodynamic  center 
of  the  NACA  0040  section  fitted  with  a  trailing  edge  wedge. 

7.4.  JfS  Section 
7.4.1.  Boundary  Layer 

Results  of  the  boundary  layer  survey  were  similar  for  both 
sections.  Figs.  7.26  and  7.27  show  the  boundary  layer  behavior  as  a 
function  of  Reynolds  number.  Transition  occurred  aft  of  the  minimum 
pressure  point  and  aft  of  the  maximum  thickness  location  at  all 
Reynolds  numbers.  Transition  occurred  at  the  30  percent  chord 
position  for  the  low  Reynolds  numbers  and  moved  toward  the  leading 
edge  to  about  the  20  percent  chord  position  at  the  higher  Reynolds 
numbers.  Separation  was  visible  only  at  Reynolds  numbers  less  than 
0.12  x  106  on  the  JfS  61M  section  and  less  than  0.165  x  106  on  the 
JfS  62M  section.  The  separation  was  of  the  laminar  boundary  layer  and 
occurred  at  the  30  percent  chord  position.  There  was  no  sign  that  it 
reattached  downstream.  The  boundary  layer  quickly  became  fully 
attached  with  a  slight  increase  in  Reynolds  number  and  remained  fully 
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7.27.  Results  of  the  boundary  layer  survey  of  the  JfS 
62M-TR-40  section. 


attached  at  all  higher  Reynolds  numbers. 

7.4.2.  Drag  Coefficient 

Fig.  7.28  shows  a  plot  of  the  measured  blockage  correction 
factors  for  the  two  JfS  modified  40  percent  sections.  The  tunnel 
blockage  ranges  from  48  percent  at  a  Reynolds  number  of  0.168  x  106  to 
12  percent  at  a  Reynolds  number  of  0.76  x  106.  Blockage  is  lower  for 
the  JfS  62M-TR-40  section  due  to  its  thinner  trailing  edge. 

The  measured  drag  coefficients,  shown  in  Fig.  7.29,  show  a  marked 
decrease  as  the  boundary  layer  became  attached  at  Reynolds  numbers 
between  0.1  and  0.2  x  106.  Above  a  Reynolds  number  of  0.4  x  10  the 
drag  coefficients  remained  approximately  constant.  As  expected,  the 
drag  coefficients  were  higher  for  the  JfS  61M  section  with  the  thicker 
trailing  edge  than  they  were  for  the  JfS  62M  section.  The  drag 
coefficient  approaches  an  approximately  constant  value  of  0.113  for 
the  JfS  61M  section  and  0.088  for  the  JfS  62M. 

Fig.  7.30  compares  the  drag  coefficients  of  the  JfS  61M-TR-40, 
JfS  62M-TR-40,  JfS  61-TR-25  and  JfS  62-TR-25  sections.  The  drag 
coefficients  for  the  JfS  25  percent  t/c  sections  were  measured  by 
Thieme  [6],  with  models  measuring  10  cm  thick  x  40  cm  chord  x  40  cm 
span  (3.9"  x  15.8"  x  15.8"),  at  the  JfS  Wind  Tunnel  at  the  University 
of  Hamburg,  Germany.  The  drag  coefficients  of  the  25  percent  sections 
were  only  slightly  less  than  those  of  the  40  percent  modified 
sections.  That  there  was  such  little  difference  indicates  the  penalty 
paid  for  the  blunt  leading  edge  of  the  unmodified  25  percent  section 
as  opposed  to  the  circular  leading  edge  of  the  modified  40  percent 
sections. 

The  base  drag  results  are  plotted  in  Fig.  7.31  for  the  two 
modified  JfS  sections.  The  fraction  of  the  total  drag  due  to  the 
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7.28.  Measured  blockage  of  the  JfS  40%  sections 
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blunt  edge  ranged  from  2  percent  to  8  percent  for  the  JfS  62M  section. 
The  ratio  of  the  area  of  the  blunt  trailing  edge  to  the  projected 
frontal  area  is  0.06.  For  the  JfS  61M  section  the  ratio  of  trailing 
edge  drag  to  total  drag  ranged  from  27  percent  to  37  percent.  The 
ratio  of  the  area  of  the  trailing  edge  to  the  projected  frontal  area 
for  this  section  is  0.20. 

For  both  of  the  modified  JfS  sections,  the  fraction  of  the  total 
drag  due  to  the  blunt  trailing  edge  was  approximately  constant  with 
Reynolds  number.  The  average  fraction  was  30.5  percent  for  the 
JfS  61M  and  15  percent  for  the  JfS  62M.  At  all  Reynolds  numbers  for 
which  trailing  edge  static  pressures  were  measured,  the  boundary 
layers  were  fully  attached.  The  measured  static  pressures  were 
negative,  i.e.  less  than  the  freestream  static  pressure. 
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7.4.3.  Hydrodynamic  Center 

The  two  JfS  modified  sections  both  show  good  weathervane 

stability.  Figs.  7.32,  7.33,  7.34  and  7.35  show  the  results  of  the 

stability  tests.  The  hydrodynamic  center  of  both  fairings  was  located 
at  the  30  percent  chord  position,  which  is  well  aft  of  the  best 
mechanical  center  located  at  the  20  percent  chord  position.  The  long, 
thick  trailing  edge  of  these  two  sections  has  the  desired  effect  of 
moving  the  hydrodynamic  aft  from  the  1/4  chord  position  predicted  by 

thin  airfoil  theory.  With  the  pivot  located  at  the  35  percent  chord 

position,  both  fairings  were  dynamically  unstable. 

Figs.  7.36  and  7.37  are  plots  of  the  moment  coefficients  for 

small  angles  of  attack  vs.  pivot  location.  This  method  also 

confirmed  that  the  hydrodynamic  center  for  both  sections  was  just  aft 

of  the  30  percent  chord  position  but  forward  of  the  35  percent  chord 

position,  approximately  0.34  chord  for  the  JfS  61  and  0.33  chord  for 

the  JfS  62  section. 
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7.32.  Yaw  equilibrium  positions  of  the  JfS  61M-TR-40  s< 
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equilibrium  positions  of  the  JfS  61M-TR-40  section 
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7.34.  Yaw  equilibrium  positions  of  the  JfS  62M-TR-40  se 
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Fig.  7.35.  Yaw  equilibrium  positions  of  the  JfS  62M-TR-40  section 
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8.0.  STUDY  CONCLUSIONS 

Each  of  the  five  sections  tested  has  hydrodynamic  characteristics 
that  recommend  its  use  as  a  fairing  for  marine  applications.  However 
none  of  the  sections  can  be  considered  ideal.  As  in  all  engineering 
problems,,  an  acceptable  solution  arises  from  a  compromise  of  the 
various  criteria.  This  section  will  compare  the  experimental  results 
of  the  five  fairings  and  evaluate  their  performance  in  terms  of  the 
desired  characteristics.  Through  this  proceedure  it  is  possible  to 
recommend  one  of  the  fairings  as  most  appropriate  for  further  study. 

8.1  Boundary  Layer 

Figures.  8.1  and  8.2  compare  the  boundary  layer  visualization 
results  for  all  five  sections  at  two  Reynolds  numbers.  The  boundary 
layer  criteria  for  a  successful  fairing  demands  a  large  region  of 
laminar  flow  extending  from  the  leading  edge  and  a  separation  free 
pressure  recovery  region  aft  of  the  maximum  thickness  point.  The 
Liebeck  fairing  was  designed  specifically  with  these  goals  in  mind. 
Of  the  five  fairings  the  Liebeck  section  had  the  largest  region  of 
laminar  flow.  The  laminar  region  was  present  for  the  leading  30  to 
47  percent  of  the  fairing  depending  on  the  Reynolds  number.  As  the 
leading  edge  profiles  increased  in  bluntness  the  extent  of  the  laminar 
boundary  layer  decreased.  The  laminar  boundary  layer  of  the  NACA  0040 
section  covered  30  to  41  percent  of  the  fairing.  Both  of  the  JfS 
modified  sections  and  the  Fathom  section  have  circular  leading  edges. 
The  laminar  to  turbulent  boundary  layer  transition  occurred  between 
the  14  and  19  percent  chord  positions  for  the  Fathom  section  and 
between  the  22  and  33  percent  chord  postions  for  the  JfS  sections. 
The  Fathom  section  had  the  smallest  extent  of  laminar  flow. 

Two  approaches  to  the  design  of  the  pressure  recovery  profile  of 
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the  sections  are  incorporated.  The  NACA  and  the  Fathom  sections  are 
convex  in  profile  aft  of  their  maximum  thickness  points,  while  the 
Liebeck  section  and  both  JfS  modified  sections  are  concave. 
Separation  occurred  at  all  tested  Reynolds  numbers  on  both  convex 
profiles  and  varied  little  with  changing  Reynolds  number.  For  all 
three  of  the  concave  profiles,  massive  turbulent  boundary  layer 
separation  occurred  at  low  Reynolds  numbers  but  became  fully  attached 
as  the  Reynolds  numbers  were  increased.  The  Liebeck  section  with  the 
largest  pressure  recovery  gradients  stayed  separated  up  to  higher 
Reynolds  numbers.  The  flow  became  fully  attached  at  Reynolds  numbers 
greater  than  0.22  x  10 6  for  the  JfS  62M,  greater  than  .25  x  106  for 
the  JfS  61M,  and  greater  than  0.40  x  10  for  the  Liebeck  section.  The 
addition  of  the  vortex  generators  improved  the  separation  resistance 
of  the  Liebeck  section.  For  this  configuration  the  boundary  layer 
became  fully  attached  at  Reynolds  numbers  greater  than  0.17  x  106, 

8.2.  Drag  Coefficient 

The  behavior  of  the  boundary  layer  affects  both  the  drag 
coefficients  and  the  dynamic  stability  of  the  fairings.  To  satisfy 
the  criteria  for  low  drag  coefficients,  the  fairing  should  have  a 
maximum  extent  of  laminar  flow  and  be  completely  free  from  separation. 
An  analysis  of  the  measured  drag  coefficients  indicates  the  importance 
of  the  boundary  layer  behavior.  Figures.  8.3  and  8.4  show  the 
measured  drag  coefficients  of  all  five  tested  fairings  normalized  with 
respect  to  projected  area  at  two  different  Reynolds  numbers. 

At  low  Reynolds  numbers  of  less  than  0.25  x  106,  where  the 


boundary  layer  was  separated  for  the  Liebeck  and  the  two  JfS  sections, 
their  drag  coefficients  were  greater  that  those  of  the  NACA  and  Fathom 
sections.  The  drag  coefficients  of  the  Liebeck  section  were  the 


highest.  The  drag  coefficients  of  the  NACA  section  were  the  lowest 
due  to  its  larger  region  of  laminar  boundary  layer  and  small  region  of 
turbulent  separation. 

For  Reynolds  numbers  greater  than  0.5  x  106 ,  when  the  boundary 
layers  were  completely  attached  for  the  Liebeck  and  jfS  sections,  the 
effects  of  the  laminar  region  and  the  blunt  trailing  edge  drag  became 
important.  The  success  of  the  Liebeck  design  is  apparent.  The  drag 
coefficients  for  this  section  were  less  that  those  of  all  of  the  other 
shapes.  Drag  coefficients  in  this  Reynolds  number  range  were  slightly 
higher  for  the  NACA  0040  shape  due  to  the  small  region  of  boundary 
layer  separation,  while  the  drag  coefficients  for  the  Fathom  section 
were  even  higher.  The  small  region  of  laminar  flow  and  the  additional 
drag  effects  of  the  blunt  trailing  edge  were  responsible.  The  two  JfS 
sections  had  the  highest  drag  coefficients  in  this  Reynolds  number 
range.  Although  they  had  a  larger  region  of  laminar  flow  then  the 
Fathom  section,  their  thickness  ratios  of  40  percent  and  their  blunt 
trailing  edges  contributed  to  their  higher  drag.  The  JfS  61M  section 
with  its  thicker  trailing  edge  had  a  much  higher  drag  than  the  JfS  62M 
section. 

It  is  important  to  note  the  high  drag  coefficients  of  the  Fathom 
section  throughout  the  tested  Reynolds  number  range.  The  Fathom 
section  has  a  t/c  of  .25  while  for  all  the  others  t/c  =  .40.  The  high 
drag  coefficients  of  this  shape  were  a  result  of  the  blunt  (circular) 
leading  edge,  the  region  of  boundary  layer  separation,  and  the 
existance  of  the  blunt  trailing  edge.  Although  the  JfS  fairings  also 
have  circular  leading  edge  profiles  and  blunt  trailing  edges,  their 
separation  free  pressure  recovery  sections  helped  to  reduce  their 
total  drag. 


Clearly  the  Liebeck  and  the  NACA  sections  best  meet  the  criteria 
of  low  drag  at  the  higher  Reynolds  numbers.  However,  the  other 
sections  have  significantly  lower  drag  coefficients  than  an  unfaired 
circular  cylinder.  At  a  Reynolds  number  of  0.6  x  106  the  drag 
coefficient  of  a  circular  cylinder  is  0.312  and  the  coefficients  based 
on  frontal  area  for  the  Liebeck,  NACA,  Fathom,  JfS  62M,  and  jfS  61M 
were  0.075,0.08,  0.136,  0.0875,  0.1325  respectively. 

8.3.  Hydrodynamic  Center 

The  most  important  criterion  for  a  successful  fairing  for  line 
structures  is  its  streamlining  or  "weathervane"  behavior.  In  order  to 
prevent  tow-off  or  "kiting"  of  the  structure  due  to  hydrodynamic  lift 
forces,  the  fairing  must  align  itself  with  the  direction  of  the  fluid 
flow.  This  criterion  requires  the  hydrodynamic  center  to  be  aft  of 
the  mechanical  center  of  rotation.  To  overcome  the  friction  forces 
between  the  fairing  and  the  line  structure  the  restoring  moments  about 
the  rotational  center  must  be  high.  Figure  8.5  depicts  these 
properties  for  the  five  sections. 

This  criterion  immediately  eliminates  the  NACA  0040  section  from 
consideration  as  an  acceptable  freely  rotating  fairing.  The  location 
of  the  hydrodynamic  center  was  found  to  be  forward  of  the  leading 
edge.  Thus,  for  all  possible  locations  of  the  rotational  center  the 
fairing  will  misalign  with  the  flow.  The  resultant  lift  forces  will 
cause  severe  tow-off  problems.  The  existance  of  a  region  of  separated 
flow  near  the  trailing  edge  at  all  Reynolds  numbers  will  also 
contribute  to  the  vibration  and  instability  of  this  fairing. 

The  Liebeck  section  will  have  similiar  problems.  The 
hydrodynamic  center  for  this  fairing  was  found  to  be  located  just 
forward  of  the  1/4  chord.  The  maximum  diameter  circular  cylinder  that 


the  Liebeck  section  can  fair  would  be  centered  at  the  35  percent  chord 
position.  For  this  configuration  the  hydrodynamic  center  would  be 
forward  of  the  mechanical  center  resulting  in  flow  misalignment  and 
kiting. 

There  are  two  possible  configurations  of  the  Liebeck  fairing  that 
could  be  used  to  overcome  this  problem.  First,  the  section  could  be 
used  to  fair  line  structures  of  smaller  diameter.  For  acceptable 
streamlining  behavior  the  maximum  diameter  must  be  less  then 
34  percent  of  the  chord  length  with  its  center  located  forward  of  the 
20  percent  chord  position.  This  is  an  inefficient  use  of  the  fairing 
shape. 

The  second  configuration  requires  the  addition  of  the  trailing 
edge  wedge.  The  hydrodynamic  center  was  now  found  to  be  located  at 
the  35  percent  chord  position,  coincident  with  the  center  of  the 
maximum  incribed  diameter.  Line  structures  with  smaller  diameters,  35 
to  39  percent  of  the  chord  length  could  be  successfully  faired.  The 
improved  streamlining  performance  is  compromised  by  the  higher  drag 
that  results  from  the  presence  of  the  wedge. 

Problems  with  kiting  of  the  prototype  Fathom  Flexnose  25  percent 
fairing  were  described  by  Henderson  [2] .  The  location  of  the 
hydrodynamic  center  was  the  cause  the  problem.  This  investigation 
confirmed  the  location  of  the  hydrodynamic  center  as  being  between 
.15c  and  ,25c  depending  on  Reynolds  number.  The  rotational  center  of 
the  Fathom  fairing  is  the  12.5  percent  chord  position.  The  measured 
restoring  moment  gradient  about  the  10  and  15  percent  positions  were 
-0.428  and  -0.215  respectively.  Henderson  measured  dC/dCT about  the 
12.5  percent  chord  position  and  found  it  to  be  -0.115.  The  low 
restoring  moment  about  the  mechanical  center  was  unable  to  overcome 
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the  frictional  moments  encountered  between  the  fairing  and  cable. 
Thus  the  fairing  would  not  maintain  alignment  with  the  flow  and  lift 
forces  would  result. 

For  both  of  JfS  sections, the  hydrodynamic  centers  were  located 
significantly  aft  of  the  best  rotational  center.  Both  sections  with 
circular  leading  edges  had  their  best  mechancical  center  located  at 
the  Xp/c  =  0.2.  The  hydrodynamic  center  for  both  fairings  was  between 
the  30  and  35  percent  chord  position.  The  restoring  moment  slopes 
about  this  mechanical  center  were  relatively  high.  The  gradient 
was  -1.09  about  xp/c  =  0.2  for  the  thin  trailing  edge  JfS  62M  section 
and  -.583  for  the  JfS  61M  section. 

With  regard  to  the  streamlining  behavior  of  the  five  fairings,  it 
was  clear  that  the  two  JfS  sections  are  superior.  The  hydrodynamic 
center  was  well  aft  of  the  rotational  center  and  their  restoring 
moment  gradients  are  high. 

8.4.  Summary 

In  summarizing  the  hydrodynmaic  characteristics  of  the  five 
fairings  it  is  evident  that  the  JfS  62M-TR-40  section  best  meets  the 
stated  criteria.  The  boundary  layer  for  this  section  was  well 
behaved.  Separation  was  a  problem  only  at  Reynolds  numbers  less  than 
.25  x  105.  Its  hydrodynamic  center  was  well  aft  of  the  ideal 
rotational  center.  The  restoring  moment  gradient  was  high.  Tow-off 
should  not  be  a  problem.  The  drag  coefficients,  although  higher  than 
those  of  the  NACA  and  Liebeck  section,  were  within  the  range  predicted 
by  Hoerner  [10]  for  other  sections  of  comparable  thickness.  The 
JfS  62M  section  with  its  thin  trailing  edge  has  lower  drag 
coefficients  then  the  JfS  61M  section.  The  JfS  61M  section  was  in 
every  other  way  comparable  to  the  JfS  62M. 
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9.0.  RECOMMENDATIONS 

The  first  phase  of  the  fairing  study  determined  that  the  JfS 
62M-TR-40  modified  section  is  worthy  of  further  study.  The  section 
utilizes  a  semi-circular  forebody  to  move  the  mechanical  rotation 
center  as  far  forward  as  possible,  a  concave  aft  section  for  boundary 
layer  separation  control  and  a  thick  trailing  edge  for  increased 
restoring  moment. 

It  is  recommended  that  the  second  phase  of  tests  continue  the 
two-dimensional  characterization  of  the  JfS  modified  section  in  the 
Venturi  tunnel  by  measuring  the  lift,  drag  and  pitching  moments  of  a 
two-dimensional  model  mounted  horizontally  between  two  vertical  plates 
to  assure  two-dimensional  flow.  The  three-component  wind  tunnel 
balance  would  be  used  to  determine  the  two-dimensional  lift  and 
pitching  moment  characteristics,  which  were  not  determined  in  the 
phase  one  tests,  and  to  determine  the  drag  coefficient  as  measured  by 
the  balance  which  would  be  used  to  check  the  wake  measurement 
technique. 

The  purpose  of  the  third  phase  of  tests,  which  would  be  conducted 
in  the  University  of  Washington  Kirsten  Wind  Tunnel,  would  be  to 
examine  the  three-dimensional  characteristics  of  the  JfS  section.  The 
Kirsten  Wind  Tunnel  is  of  a  closed  circuit,  double  return  type,  with  a 
2.44  by  3.66  m  test  section  vented  to  the  atmosphere.  Wind  velocities 
up  to  111/75  m/sec  can  be  generated  in  the  test  section,  corresponding 
to  dynamic  pressures  of  1  to  160  psf.  A  30.5  cm  model  chord  would 
allow  tests  in  the  Kirsten  wind  tunnel  to  a  maximum  Reynolds  number  of 
2.5  x  106.  A  vertically  mounted  model  comprising  several  discrete 
segments  of  the  fairing  section  elements,  Fig.  9.1,  would  be  mounted 
in  the  tunnel.  The  30.5  cm  model  would  allow  approximately  eight 


segments  of  the  fairing  to  be  simulated,  assuming  an  aspect  ratio  of 
one  for  each  fairing  segment.  It  should  be  noted  that  aspect  ratios 
of  two  would  also  be  investigated.  An  instrumented  floating  fairing 
element,  located  at  approximately  midspan  of  the  model,  would  be  used 
to  measure  the  following: 

(1)  The  element  section  would  be  built  with  a  distribution  of 
chordwise  pressure  taps  to  measure  the  chordwise  pressure 
distribution. 

(2)  The  upper  and  lower  ends  of  the  section  would  be  mounted  with 
block  gauges.  Fig.  9.1,  two  of  which  would  measure  the  normal 
force  component  and  a  third  mounted  at  the  top  which  would 
measure  the  tangential  forces. 
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DATA  REDUCTION  PROGRAM 
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DRAG  COEFFICIENT 
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REAL  M,  INTEQ( 25) 

DIMENSION  Y(25),Q(25),QRATIO(25) 

BYTE  PARAMf 11), INFIL( 11), DATFIL{ 1 1 ) 

BYTE  XCONF ( 40 ) , XNAME ( 6 ) , XDATE( 20 ) ,NRUN ( 2 ) 

FORMAT( 2X , *  TO  RUN  PROGRAM  USING  OLD  DATA  TYPE  -  ”0'"  ,/, 

2X, ' TO  COLLECT  NEW  DATA  TYPE  -  "N"',/, 

2X , 1 ( O/N ) ?  -') 

FORMAT (2X, 'TYPE  FILENAME  CONTAINING  PARAMETER  DATA  -') 

FORMAT(A2) 

FORMAT(2X, 'INPUT  RUN  NUMBER  -  ') 

FORMAT(2X,' INPUT  QI,  HIGH  OR  LOW  SCALE',/, 

2X, 'FORMAT(F6.0,H/L)-' ) 

FORMAT(2X, 'INPUT  ROOM  TEMPERATURE  (C)  -') 

FORMAT(F6.0,A2) 

FORMAT(F6.0) 

FORMAT( 2X , ' OPEN  VALVE  NUMBER:  ' , 13 , / , 2X , ' TYPE  Y  TO  CONFIRM  -  ') 
FORMAT( 2X , ' TAP  NO . ' , 5X , ' Y ( IN . ) ' , 12X,'Q(PSF) ' ,/) 

FORMAT( 4X, 13 ,8X, F5. 2 , 8X , F 1 0 . 4 ) 

FORMAT(2X, 'TYPE  THE  TWO  TAP  NUMBERS  THAT  DEFINE  THE  EDGE  ',/, 

2X , ' OF  THE  WAKE.  NEGATIVE  TAP  #  FIRST:  ') 

FORMAT(//,2X, 'DO  YOU  WISH  TO  CONTINUE  DATA  REDUCTION  [Y/N]?  -') 
FORMAT( 213 ) 

FORMAT(2X, 'TYPE  FILENAME  OF  OLD  DATA  -*) 

FORMAT(2X, 'TYPE  FILENAME  FOR  DATA  STORAGE  -’ ) 

FORMAT(2X, 'WOULD  YOU  LIKE  A  PRINTED  COPY  OF  THE  RESULTS  [Y/N]? 

-  ') 

FORMAT( 2A1 ) 

FORMAT ( 1 1A1 ) 

FORMAT( 2X, ' Q( I )* ' ,F10.5, *  (PSF)',/) 

FORMAT(2X,  'OPEN  VALVE  FOR  T.E.  STATIC  PRESSURE  TAP',/, 

2X , ' TYPE  Y  TO  CONFIRM',/, 

2X , ' TYPE  N  IF  NO  TE  PRESSURE  TAP  [Y/N]  -  ') 

FORMATf  2X, ' Q  STATIC  AT  T.E.-' ,F10.5, '  (PSF)',/) 

FORMATf  2X, 'DO  YOU  WANT  TO  INPUT  THE  SEPARATION  LOCATION  [Y/N]?  • 

) 

FORMATf 2X, 'INPUT:  XSEP/C  -  ') 

FORMAT( F5. 3 ) 

FORMAT( 2X , ' DO  YOU  WISH  TO  CHANGE  TRANSDUCER  CALIBRATION  [Y/N]  ') 
FORMATf  2X, ' INPUT  "Q“M*VBAR+B" ;  M,B-  ') 
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m 

6000 

FORMATS  6A1 , 20A1 ) 

w 

6001 

FORMAT( 6A 1 , 2  0A 1 , 3F 1 5 . 8 ) 

6101 

FORMAT( 11A1,20A1,3F15.8) 

6005 

FORMAT{ 40 A 1 ) 

ft 

6010 

FORMAT* 3F 1 0 • 4 ) 

K. 

6011 

FORMAT* 3F7.4) 

6019 

FORMAT( 2F 10. 4, 12 ) 

w 

6020 

FORMAT{ 2F10.4,I2,F5.3) 

6030 

FORMAT( 13) 

6040 

FORMAT( F10 . 5 ) 

r*. 

6050 

FORMAT{ 3F 1 5 . 10) 

c: 

7000 

FORMAT* 13) 

7010 

FORMAT* 2F1 0.5) 

9000 

FORMAT* 2X, 'END  OF  PROGRAM',/) 

S 

9001 

FORMAT* 2X, 'DO  YOU  WISH  TO  RUN  THE  PROGRAM  AGAIN 

C 

C 

V. 

lJI 

C 

PART  1 

£ 

C 

DETERMINE  IF  OLD  DATA  IS  TO  BE  RECALCULATED 

C 

OR  IF  NEW  DATA  IS  TO  BE  TAKEN 

£ 

C 

m=o.o 

B=0 . 0 

« 

1 

CALL  HOME 

WRITE* 5, 5000) 

READ*  5,5010)  IC 

IF(IC.NE. 'N' )  GO  TO  500 


READ  FILE  CONTAINING  OPERATING  PARAMETERS 


WRITE( 5, 5005) 

READ( 5,5075) PARAM 
CALL  OPEN ( 6 , PARAM, 2 ) 

READ( 6,6000 )  XNAME , XDATE 

READ (6, 6005)  XCONF 

READ (6, 60 10 )  F,C.T 

READ ( 6 , 60 1 1 ) XTW , TWT, TWW 

READ (6, 60 19)  THETA, XRAXE,NTAPS 

READ(6,6030)  NPTS 

DO  10  I«1,NTAPS 

INTEQ(I)«0.0 

READ* 6, 6040)  Y(I) 

ENDFILE  6 


GET  PARAMETERS  FROM  TERMINAL 


WRITE ( 5,5019) 

READ( 5,5074 )NRUN 
WRITE( 5,5020 ) 

READ( 5, 5030)  QI,IS 


>A.SV>V 


,  a  Ji_*  w 7  6  .m 


kl  |U  . 


WRITE( 5,5025 ) 

READ( 5,5033 )  TC 

RHO*. 002522/ ( 1 .0+.00367*TC) 

R MU*. 14248E-8*TC+. 35242E-6 
QTE=0 . 0 

CALCULATE  QU  -  UPSTREAM  DYNAMIC  PRESSURE 

IF{ IS. EQ. ' H' )  GO  TO  30 

QU-.96*QI 

GO  TO  40 

QU= . 983*QI- . 24 

CONTINUE 

DATA  COLLECTION 

WRITE( 5,5100) 

READ( 5, 5074 ) IB 
IF(IB.NE. ’Y' )  GO  TO  49 
WRITE(5,51 10) 

READ( 5,7010)M, BS 
IF( M.EQ. 0 . 0 )GO  TO  48 
NN=1+(NTAPS-1)/2 
DO  59  I“1 ,NTAPS 
J-I-NN 

WRITE(5,5035)  J 

READ( 5,5010)  IT 

IF( IT.NE . 'Y' )  GO  TO  50 

DO  54  KK=1 , 500 

DUM=SQRT(  (BS*M)*(BS*M)  )/2. 

CONTINUE 

CALL  ADCONV( NPTS, VBAR) 

Q(I)«M*VBAR+BS 
WRITE( 5, 5080 )  Q(I) 

CALL  BEEP( 50 , 200 ) 

IF(Q(I) .LT.0.0)  GO  TO  50 

CONTINUE 

WRITE(  5 , 5090 ) 

READ( 5,5010) IT 
IF( IT.EQ. ' Y' )GO  TO  65 
IF( IT. EQ. *N' )GO  TO  70 
GO  TO  60 

CALL  ADCONV(NPTS, TEV) 

QTE=  M*  TEV+B  S 
WRITE ( 5,5091 )QTE 
WRITE(5,5095) 

READ (5,5010) IX 
X SEP-0.0 

IF( IX. NE. ' Y' )  GO  TO  100 
WRITEf  5, 5096 ) 


8 

I 

8 


s 


& 


i 


& 


<* 


.<• 


* 


,*  • 

A 
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READ( 5 , 5098 )X SEP 
100  CONTINUE 

C 

C  PLOT  AND  LIST  DATA  ON  TERMINAL 

C 

CALL  SUBPLT(NTAPS,Y,Q) 

WRITE( 5, 5040 ) 

DO  110  1-1, NTAPS 
J-I-NN 

110  WRITE<5,5050)  J,Y(I),Q(I) 

WRITE { 5 , 5060 ) 

READ(5,5010)  II 

IF(I1 .EQ. 'N' )  GO  TO  9990 

WRITE ( 5,5055) 

READ(5,5065)  K1,K2 

K1-K1+NN 

K2-K2+NN 

CO“(C(K1)+Q(K2))/2.0 

C 

C  SAVE  DATA  ON  FILE 

C 

IF( IC.NE. 'N' )  GO  TO  600 
DO  115  1-1,6 

115  DATFIL(I)-XNAME(I) 

DATFIL(7)-NRUN< 1 ) 

DATFIL( 8 ) — NRUN { 2 ) 

DATFIL(9)-'D' 

DATFIL( 1 0 )  —  *  A ' 

DATFIL( 1 1 )-'T' 

CALL  OPEN ( 8 , DATF IL , 2 ) 

WRITE(8,6101 )DATFIL,XDATE,QU,QI,QTE 
WRITE( 8 , 6005 ) XCONF 
WRITE(8,6010)B,C,T 
WRITE(  8,6011 )XTW,TWT,TWW 
WRI TE( 8, 6020 ) THETA, X RAKE, N TAPS, XSEP 
WRITE( 8,6030 )NPTS 
WRITE(8,6050)TC, RHO,  RMU 
WRITE(  8 , 7000  )  NTAPS 
DO  120  1-1, NTAPS 
120  WRITE( 8,7010)  Y(I),Q(I) 

ENDFILE  8 
GO  TO  600 
C 

C  GET  OLD  DATA  FROM  FILE  IF  REQUESTED 

C 

500  WRITE(  5, 5070  ) 

READ( 5 , 5075 )  INFIL 
CALL  OPEN( 7, INFIL, 2 ) 

READ ( 7,6101 )DATFIL,XDATE,QU,QI,QTE 
READf  7, 6005) XCONF 


100 
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j 
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READ {7,6010)B,C,T 

READ ( 7 , 60 1 1 ) XTV , TWT , TWW 

READ ( 7 , 6 0 2 0 ) THETA , X  RAKE , NTAP S , X SEP 

READ( 7,6030 )NPTS 

READ (  7,6050  )TC,  RHO,  RMU 

READ( 7, 7000 )  NTAPS 

DO  510  1=1, NTAPS 

INTEQ( I ) =0 . 0 

510  READ (7,7010)  Y(I),Q(I) 

ENDFILE  7 
NN=1+(NTAPS-1 )/2 
GO  TO  100 
C 

C  CALCULATION  OF  FLOW  PARAMETERS 

C 

600  VU=SQRT(2.0*QU/RHO) 

VO=  SQRT  (  2 . 0  *  QO/ RHO ) 

E= ( VO-VU ) /VU 
REN=C*  RHO* VO/ ( RMU* 12.) 

RENU=C*  RHO* VU/ ( RMU* 12.) 

CALCULATION  OF  INTEGRAND  FOR  DRAG  INTEGRATION 
C 

DO  650  I=K1,K2 
QRATIO( I )=Q( I )/QO 

650  INTEQ(I)=SQRT(QRATIO(I) )-QRATIO(I) 

CALL  INTGRL  ( Y, INTEQ, NTAPS, AREA ) 

C 

C  CALCULATE  DRAG  COEFFICIENTS 

C 

CDBC=(2./C)*AREA 

CDBT=CDBC*C/T 

C 

C  WRITE  DATA  AND  RESULTS 

C 

L=5 

CALL  HOME 
IP= ' N ' 

675  WRITE(L, 1790) 

WRITE(L, 1800 )XCONF, XDATE 
WRITE( L, 1801 )DATFIL 

WRITE(L, 1805 )QU, QI , RENU, TC, VU, T, RHO, C, RMU, B, E, XRAKE , QTE 
IF(XTW.EQ.O.O)  GO  TO  699 
WRITE( L, 1806)XTW 
WRITE{ L, 1807)TWT,TWW 
699  IF(L.EQ.5)PAUSE 

WRITE( L,  1810) 

DO  700  1=1, NTAPS 
J=I-NN 

WRITE(L, 1820)  J,Y(I),Q(I), INTEQ ( I ) 


700 


WRITE(L, 1835 ) REN 
WRITE(L, 1830)  CDBC,CDBT 
WRITE(L, 1840) 

IF { X SEP • EQ >0.0) GO  TO  720 
WRITE(L,  1842  )XSEP 
GO  TO  730 
WRITEfL, 1843) 

IF(L.EQ.2)  WRITE(2, 1831 ) 

WRITE( 5, 5072) 

READ( 5, 50 10 ) IP 

IF( IP.NE. ' Y' )  GO  TO  9990 

L-2 

GO  TO  675 

FORMAT( 5X, ' SYMMETRICAL  FAIRING  SECTIONS' ,/, 5X, 

1  '2-DIMENSIONAL  DRAG  COEFFICIENT  TESTING') 

FORMAT(/,5X, 'MODEL  CONFIGURATION  -  * , 40A1,/,5X, 'DATE  -  \20A1) 
FORMAT( /,  5X,  ' TEST  DATA  FILENAME  -',11A1) 

FORMAT(5X, 'LOCATION  OF  BOUNDARY  LAYER  TRIP  WIRE,  (XTW/C)-  ',F4.2 
1  ) 

FORMAT( 5X , ' TRIP  WIRE  DIMENSIONS:  THICKNESS  ',F7.4,'  WIDTH-  ', 

1  F7.4, '  INCHES',/) 

FORMAT(/,9X, 'TAP  NO.',7X,'Y  (IN.)',7X,'Q  (PSF) ' , 7X, 'DRAG 
1 INTEGRATION' ,/, 

1  54X,' INTEGRAND', /, 

29X,  ' - *  ) 

FORMAT( 10X,I3,2X,3(5X,F10.5) ) 

FORMAT(/,5X, 'CD(BC)-' ,F10.4,/, 5X, 'CD{BT)=' ,F10.4) 

FORMAT ( 1H1 ) 

FORMAT(//,5X, 'DOWNSTREAM  REYNOLDS  NUMBER;  REN=  *,E8.3) 
FORMAT(//,5X, 'SEPARATION  LOCATION  AS  INDICATED  BY  MICROTUFT 
1 BEHAVIOR: ' ) 

FORMAT( 5X, ' XSEP/CHORD  -  ’,F5.3,//) 

FORMAT(5X, ' SEPARATION  LOCATION  WAS  NOT  RECORDED  FOR  THIS  RUN') 
FORMAT( // , 5X, ' QU— ' , F8 .4, IX, ' (PSF) ' , 14X, 'QI-' ,F8.4, IX, ' (PSF) ' ,/, 

1  5X, 'UPSTREAM  REN=  ', E8. 3, 9X, ' TEMPERATURE- ' ,F6. 2, IX, ' (CENTIGRADE) 
1  ' 

2  5X, 'VEL-' ,F10.3, IX, ' (FT/ SEC) ’ , 8X, ' THICKNESS® ' ,F6.2, IX, ' (IN. ) ’ ,/, 

3  5X,  'RHO-'  ,F10.8,  IX,  '  (TE  UNITS)  '  ,6X, 'CHORD-' ,F6.2,  IX, '  (IN.  )  '  ,/, 

4  5X, 'MU-' ,F12. 10, IX, ’ (TE  UNITS) ' , 5X , ' SPAN-' ,F6. 2, IX, ' ( IN. ) ' ,/, 

5  5X, 'BLOCKAGE-' ,F8. 5, 14X, 'DISTANCE  TO  WAKE  RAKE-' , F6. 2, IX, ' ( IN. ) ' 

6  ,/ , 36X, ' T. E.  STATIC  PRESSURE-' ,F8. 4, '  (PSF)',/) 


9990 

WRITE( 5,9001 ) 

READ( 5, 50 1 0 ) IA 

v: 

& 

IF( IA.NE. 'N‘ )  1 

9999 

WRITE( 5,9000) 
STOP 

END 

SUBROUTINE  FOR  A/D  CONVERSION 
OF  10  VOLT  SIGNAL  ON  CHANNEL  0 
OF  VD  BOARD 

SUBROUTINE  ADCONV(NPTS, VBAR) 
INTEGER  A,A1 , A2, A3 


* 

5; 


A-57536 
Al-A+1 
A2— A+2 
A3-A+3 

CALL  POKE( A, 252 ) 

SUM-0. 

CONVERSIONS  ON  CHANNEL  0 
CALL  POKE ( A 1 , 0 ) 

START  CONVERSIONS 

DO  500  I«1,NPTS 
CALL  POKE ( A2 , 0 ) 

TEST  FOR  END  OF  CONVERSION 

IF(IPEEK( A1J.LT. 128)  GO  TO  10 
READ  REGISTERS 

L— IPEEK( A2) 

M*IPEEK(A3) 

D«256.*M*L 

IF(D*GT. 2047 . )  GO  TO  50 

PSUf*=D/NPTS 

GO  TO  500 

PSUM=( D-4095 )  /NPTS 
D  SUM=SUM+PSUM 

VBAR=( SUM* 10.) /2047 . 

WRITE ( 5 ,5001)  VBAR 

01  FORMAT( 2X , ' VBAR- ' ,  F 1 0 . 5  ) 

RETURN 

END 

SUBROUTINE  SUBPLT(NTAPS,X,Q) 

REAL  X(NTAPS) ,Q(NTAPS) 

INTEGER  XPT, YPT 

WRITE(5, 1000) 

00  FORMAT(2X, 'WELCOME  TO  SUBPLOT,  TYPE  <CR>  TO  CONTINUE',/) 

PAUSE 

NT1-NTAPS-1 
QMAX«Q( 1 ) 

DO  100  1-1, NT1 

IF(Q(I)  .LE.Q<I+1)  )  QMAX-Qd+1) 

0  CONTINUE 

YSCL-AB S( Q MAX/4 0 . ) 

IF ( QMAX . EQ .0.0)  YSCL-1. 


ft 


I 


f-  *.w 


■  »-*  Ik'  l_t 


CALL  GR( 1,0) 

CALL  COLOR( 15) 

CALL  VLIN( 0,47,0) 

CALL  HLXN(0,39,47) 

CALL  COLOR( 7 ) 

DO  300  1*1 ,NTAPS 
XPT*19+INT(X( I )/. 25 ) 

IF( XPT. LE. 0 .OR. XPT. GE.40 )  GO  TO  300 
YPT*48-INT ( Q { I ) /Y  SCL ) 

CALL  VLIN(YPT,46,XPT) 

CONTINUE 
PAUSE 
CALL  ROME 
CALL  TEXT 
RETURN 
END 


c 


SUBROUTINE  INTGRL(X, Y,N, AREA) 


DIMENSION  Z(5) ,XP( 10) ,F(5) ,XI(5) 

DIMENSION  DY ( 1 1 0 ) 

'  REAL  X(N),Y(N) 

NINT* 10 
II-1 

;  AREA-0 . 0 

Z(1)-0.0 

f  Z( 2 )— 0 . 53S4693 1 

'  Z(  3)— Z(  2 ) 

Z(4)-0. 90617985 
Z(5)— Z(4) 

W1-0. 56888889 
W2-0. 47862867 
W4-0 . 23692689 
W3-W2 
W5-W4 

CALL  MONDER(X,Y,N,DY) 

DX— ( X(N) -X( 1 ) )/NINT 
DO  10  I—  1 , NINT 

10  XP(I)-X( 1 )+FLOAT(I-1 )*DX 

DO  20  1-1, NINT 
XC-XP(I)+DX/2. 

DO  15  J-1,5 
XI(J)-XC+Z( J)*DX/2. 

CALL  PWCFEV(II,N,X,Y,DY,II,XI( J) ,F( J) ) 

XINT-(W1*F( 1)+W2*F(2)+W3*F(3)+W4*F(4)+W5*F(5) )*DX/2. 
20  AREA-AREA+XINT 

RETURN 
END 


C 

SUBROUTINE  MDNDER(X, Y,N,DY) 

C  IN  THIS  VERSION,  ONLY  25  POINTS  CAN  BE  FIT 
C 

DIMENSION  X(25),Y(25) , DY( 25) , OCH( 25) 

C 

C  CHECK  THAT  THERE  ARE  MORE  THAN  3  POINTS 
C 

IF(N.GT. 3 )  GO  TO  20 
TYPE  1100, N 

1100  FORMAT (//IX* ERROR— MONDE R  REQUIRES  AT  LEAST  4  DATA* 
1*  POINTS— N-  '14/  ) 

STOP 

20  CONTINUE 
C 

C  EVALUATE  THE  CHORD  LENGTHS 
C 

IERR-0 
OCH( 1)-0.0 
DO  10  1-2, N 
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DC-X(I)-X(I-1) 

IF(DC.LE.O.O)  IERR*1 
OCH(I)**OCH(I-1)+DC 
10  CONTINUE 

XF( IERR.EQ. 0)  GO  TO  30 
TYPE  1110 

1110  FORMAT( //IX' ERROR—  MDNDER  REQUIRES  THAT  THE  X  VALUES  * 
1* STRICTLY  INCREASE'/) 

STOP 

30  CONTINUE 


NOW  CALL  MOND1  TO  GENERATE  THE  DERRIVATIVES  AT  THE  DATA  POINTS 


C 

C- 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

o 

c 

c 

c 

c 

c 

c 

c 


CALL  KJND1  (N,  OCH,  Y,  DY) 

RETURN 

END 

SUBROUTINE  MOND1  (N,  X,  F,  D) 

INTEGER  N 

REAL  X(25),  F( 25) ,  D<25),  SLOPE(25),  H{25) 


MDNDl  USES  THE  FRITSCH-CARLSON  FORUMLAS  TO  SET  DERIVATIVE 
VALUES  FOR  A  PIECEWISE  CUBIC  INTERPOLANT  TO  THE  DATA  (X,  F) 

SO  THAT  THE  INTERPOLANT  IS  MONOTONE  ON  ANY  SUBINTERVAL  ON 
WHICH  THE  DATA  ARE  MONOTONE. 

THIS  VERSION  USES.  . 

1.  THREE-POINT  DIFFERENCE  FORMULAS  TO  INITIALIZE  DERIVATIVES 
(INCLUDING  ENDPOINTS). 

2.  REGION  S( 3) . 

3.  ALGORITHM  A  FOR  MOVING  A  POINT  INTO  REGION. 

4.  ANY  NEGATIVE  ALPHA  OR  BETA  (INDICATING  A  CHANGE  IN  MONO¬ 
TONICITY  OF  THE  DATA)  IS  SET  TO  ZERO  TO  INSURE  THE  STRICT 
PIECEWISE  MONOTONICITY  OF  THE  INTERPOLANT. 

SUBROUTINE  PWCFEV  MAY  BE  USED  TO  EVALUATE  THE  RESULTING 
PIECEWISE  CUBIC  FUNCTION. 

REFERENCE..  F.  N.  FRITSCH  AND  R.  E.  CARLSON,  PIECEWISE  CUBIC 
INTERPOLATION  METHODS,  LAWRENCE  LIVERMORE  LABORATORY  REPORT 
UCRL-81230  (NOVEMBER  1978). 


ON  INPUT. . 

N  IS  THE  NUMBER  OF  DATA  POINTS. 

RESTRICTION..  N.GE.4  (NOT  CHECKED) . 

X  IS  THE  ARRAY  OF  INDEPENDENT  VARIABLE  VALUES. 

RESTRICTION..  X  MUST  BE  STRICTLY  INCREASING,  THAT  IS 
X(I)  .LT.  X(I+1),  I-MDN-1  (NOT  CHECKED). 
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IS  THE  ARRAY  OF  DEPENDENT  VARIABLE  VALUES. 


ON  OUTPUT. . 

D  WILL  BE  SET  TO  THE  DESIRED  DERIVATIVE  VALUES. 

H  WILL  BE  THE  ARRAY  OF  INTERVAL  LENGTHS, 

H(  I )  -  X(I+1)  -  X(  I)  /  I-KDN-1. 

SLOPE  WILL  BE  THE  ARRAY  OF  SLOPES  OF  CHORDS, 

SLOPE(I)  -  (F(  1+1 )  -  F(  I )  )/H(  I ) ,  I-KDN-1. 


NOTE. .  ARRAYS  H  AND  SLOPE  ARE  NO  LONGER  NEEDED  AFTER  THE  CALL  TO 
M0ND1. 


FORTRAN  INTRINSICS  USED.  .  ABS. 


ALGORITHM  BY..  F -  N.  FRITSCH,  LAWRENCE  LIVERMORE  LABORATORY,  AND 
R.  E.  CARLSON,  GROVE  CITY  COLLEGE,  PA. 

PROGRAMMED  BY. -F.  N.  FRITSCH. 

DATE  LAST  CHANGED. .  11  JANYARY  1979  (FNF) 


CHANGE  RECORD. . 

78-12-07  MINOR  COSMETIC  CHANGES  TO  GET  READY  FOR  LIBRARY. 

78- 12-20  1.  REMOVED  ARGUMENT  ICOUST. 

2.  CHANGED  ARGUMENT  Y  TO  F  (TO  BE  CONSISTENT  WITH 
PWCFEV) . 

79- 01-11  1.  CHANGED  TREATMENT  OF  INTERVAL  ADJACENT  TO  CHANGE 

IN  MONOTONICITY  OF  DATA  (SEE  ITEM  4,  ABOVE). 

2.  MINOR  ADDITIONS  TO  COMMENT  SECTION. 


LOCAL  DECLARATIONS. 


INTEGER  I,  NLESS1 

REAL  ALPHA,  BETA,  DELTA,  FUZZ,  TAU 
DATA  FUZZ  /1.0E-14/ 


INITIALIZE. 


NLESS1  =  N  -  1 


COMPUTE  INTERVAL  LENGTHS  AND  SLOPES. 


DO  10  1-1,  NLESS1 

H( I )  -  X(I+1)  -  X( I ) 

SLOPE ( I )  -  (F( 1+1 )  -  F( I ) )/H( I ) 
10  CONTINUE 


INITIALIZE  D( 1 )  VIA  NON-CENTERED  THREE-POINT  FORMULA. 


V  V* 


r 


c 

c 

c 


D(1)  -  ((H(1)+H(1)+H(2))*SL0PE(1)  -  H( 1 )*SL0PE(2) )/(H( 1 )+H(2) ) 
IF  ( D ( 1 ) *  SLOPE ( 1 )  .LT.  0.)  D(1)  -  0. 


I 


1 

§ 


I 


% 


a 


8 

$ 


CYCLE  THROUGH  ALL  INTERVALS. 


C 

C 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 


C 

C 

C 


c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 


DO  50  1-1,  NLESS1 

IF  (I  .LT.  NLESS1)  GO  TO  20 

SPECIAL  CASE  OF  RIGHT  ENDPOINT. 

D(N)  -  ( (H(N-1 )+H(N-1 )+H(N-2) )* SLOPE (N-1 ) 

*  -  H(N-1)  * SLOPE (N-2) )/(H(N-2)+H(N-1 ) ) 

IF  ( D ( N ) *  SLOPE ( N- 1 )  .LT.  0.)  D(N)  -  0. 

GO  TO  25 
20  CONTINUE 

USE  THREE-POINT  FORMULA  TO  INITIALIZE  RIGHT-HAND 
DERIVATIVE  FOR  INTERVAL  (X(I),  X(I+1>)  . 

D(  1+1 )  -  (H(I+1 )* SLOPE ( I)  +  H( I )* SLOPE( 1+1 ) )/( H( I )+H( 1+1 ) ) 

25  CONTINUE 

ADJUST  D(I)  AND/OR  D(I+1),  IF  NECESSARY  TO  INSURE  MONOTONICITY 
ON  INTERVAL  (X(I),  X(I+1))  . 

TAKE  CARE  OF  FLAT  DATA. 

IF  ( AB S( SLOPE ( I ) )  .GT.  FUZZ)  GO  TO  30 
ALPHA  -  0. 

BETA  -  0. 

GO  TO  45 
30  CONTINUE 

COMPUTE  SCALED  DERIVATIVES. 

ALPHA  -  D( I )  /  SLOPE ( I ) 

BETA  -  D{ 1+1 )  /  SLOPE ( I ) 

TAKE  CARE  OF  NONMONOTONE  DATA. 

ASSERTION..  IF  EITHER  OF  THE  FOLLOWING  TESTS  IS  SATISFIED, 
(ALPHA, BETA)  IS  NOT  IN  FIRST  QUADRANT,  WHICH 
MEANS  THAT  SLOPE  CHANGES  SIGN  AT  ONE  OR  BOTH  ENDS 
OF  INTERVAL. 

IF  (ALPHA  .LT.  0.)  ALPHA  -  0. 

IF  (  BETA  .LT.  0.)  BETA  -  0. 

ASSERTION..  ALPHA  AND  BETA  ARE  NOW  BOTH  NONNEGATIVE. 

MAKE  ALPHA  +  BETA  .LE.  3  . 


DELTA  -  ALPHA  +  BETA 
IF  (DELTA  .LE.  3.)  GO  TO  45 

ASSERTION..  POINT  IS  OUTSIDE  THE  TRIANGLE.  NEED  TO  ADJUST. 
TAU  -  3. /DELTA 
ALPHA  -  TAU* ALPHA 
BETA  -  TAU* BETA 

RECOMPUTE  DERIVATIVE  VALUES. 

45  CONTINUE 

D(I)  *  ALPHA*  SLOPE ( I ) 

D( 1+1 )  -  BETA*  SLOPE ( I ) 

50  CONTINUE 

END  OF  DERIVATIVE  ASSIGNMENT. 

RETURN 

END 

SUBROUTINE  PWCFEV  (IDERIV,  N,  X,  F,  D,  NE,  XE,  FE) 

INTEGER  IDERIV,  N,  NE 

REAL  X(N) ,  F(N) ,  D(N),  XE(NE),  FE(NE) 


EVALUATE  THE  FIRST  (IDERIV-1)  DERIVATIVES  OF  THE  PIECEWISE 
CUBIC  FUNCTION  DEFINED  BY  N,  X,  F,  D  AT  THE  POINTS  XE(I), 
I  -  1(1)NE. 


ON  INPUT. . 
IDERIV 


ON  RETURN. 


INDICATES  HOW  MANY  DERIVATIVES  ARE  DESIRED. 
RESTRICTION.-  1  -LE.  IDERIV  .LE.  3  (NOT  CHECKED) . 

NOTE..  IDERIV-1  IMPLIES  ONLY  FUNCTION  VALUES  REQUESTED- 
IS  THE  NUMBER  OF  DATA  POINTS. 

RESTRICTION.-  N  .GE.  2  (NOT  CHECKED). 

IS  THE  ARRAY  OF  INDEPENDENT  VARIABLE  VALUES. 

THE  SEARCH  PROCEDURE  ASSUMES  THAT  X  IS  STRICTLY 
INCREASING.  (NOT  CHECKED) 

IS  THE  CORRESPONDING  ARRAY  OF  FUNCTION  VALUES. 

IS  THE  CORRESPONDING  ARRAY  OF  DERIVATIVE  VALUES. 

IS  THE  NUMBER  OF  POINTS  AT  WHICH  EVALUATION  IS  DESIRED. 
IS  THE  ARRAY  OF  EVALUATION  POINTS. 

THE  SEARCH  PROCEDURE  ASSUMES  THAT  XE  +S  MDNOTONE 
INCREASING.  (NOT  CHECKED) 


CONTAINS  THE  FUNCTION  VALUES,  AS  FOLLOWS. 

FE( I , J)  IS  THE  VALUE  OF  THE  (J-I)-ST  DERIVATIVE  OF 
THE  PIECEWISE  CUBIC  AT  XE(I),  I«1(1)NE, 


LOCAL  DECLARATIONS. 

INTEGER  I,  IER,  IL 
REAL  DX,  H(4),  RDX,  T 

MAIN  EVALUATION  LOOP. 

IL  -  0 

DO  50  I  -  1,  NE 

LOCATE  INTERVAL  CONTAINING  XE(I). 

CALL  SEARCH  (N,  X,  XE(I),  IL,  IER) 

EVALUATE  HERMITE  BASIS  FUNCTIONS  AND  NEEDED  DERIVATIVES. 

CALL  HBASEV  ( IDERIV,  X(IL),  X(IL+1),  XE(I),  H) 

EVALUATE  CUBIC  AND  APPROPRIATE  DERIVATIVES. 

DO  40  ID  *  1,  IDERIV 

FE( I )  -  F( IL)*H( 1 )  +  F( IL+1 )*H( 2 ) 

*  +  D( IL) *H( 3 )  +  D(IL+1)*H(4) 

40  CONTINUE 

50  CONTINUE 

RETURN 

END 

SUBROUTINE  SEARCH  (N,  X,  XVAL,  IL,  IER) 

INTEGER  N,  IL,  IER 
REAL  X(N) ,  XVAL 


SEARCH  FOR  XVAL  IN  ARRAY  X. 

PERFORMS  A  LINEAR  SEARCH,  FROM  LEFT  TO  RIGHT.  FOR  IMPROVED 
EFFICIENCY  WHEN  LOCATING  AN  INCREASING  SEQUENCE  OF  XVAL" S, 

THE  STARTING  INDEX  FOR  THE  SEARCH  MAY  BE  SPECIFIED  BY  THE  USER. 


ON  INPUT. • 

N  IS  THE  NUMBER  OF  POINTS  IN  ARRAY  X. 


.  •  •  '  V  •/  *.*  %'  '/  s’ 


■ vV  -\v.v 


PROGRAM  TORQUE 


DATA  REDUCTION  PROGRAM 
FOR  THE  CALCULATION  OF  THE 
TORQUE  COEFFICIENT 
FROM  WIND  TUNNEL  DATA 

WRITTEN  BY: 

DAVID  GRAY 
FEBURARY  1982 

DIMENSION  Y( 25 ) , TM( 25 ) , CM( 25 ) 

REAL  M 

BYTE  XCONF(40) ,XNAME( 6) ,XDATE(20) 

BYTE  PARAM( 1 1 ) 

FORMAT(2X, 'RESTORING  MOMENT  DATA  COLLECTION  PROGRAM',/, 

2X, 'BE  CERTAIN  THAT  THE  ANALOG  SIGNAL  TO  A/D  IS  INPUT  TO 
CHANNEL  6* ,/,2X, 'TYPE  <CR>  TO  CONTINUE  -  ') 

FORMAT(/,2X, 'TYPE  FILENAME  OF  PARAMETER  DATA  -  ') 

FORMAT(A2) 

FORMAT(2X, ' INPUT  RUN  NUMBER  -  ') 

FORMAT(  13 ) 

FORMAT(2X, 'INPUT  QI,  HIGH  OR  LOW  SCALE',/, 

2X, ' FORMAT (F6.0,H/L)-' ) 

FORMAT(2X, 'INPUT  ROOM  TEMPERATURE  (C) :  ') 

FORMAT( F6 • 0 , A2 ) 

FORMAT(  F6. 0 ) 

FORMAT(2X, 'WOULD  YOU  LIKE  A  PRINTED  COPY  OF  THE  RESULTS  [Y/N] ? 

-  ') 

FORMAT(  1 1A1 ) 

FORMAT( 2X , ' DO  YOU  WISH  TO  CHANGE  TRANSDUCER  CALIBRATION  [Y/N]  ') 

FORMAT ( 2X , ' INPUT  FORCE- M* VBAR+B f  M,B-  ’) 

FORMAT( 2F 1 0 . 5 ) 

FORMAT( 6A 1 , 2  OA 1 ) 

FORMAT( 40A1 ) 

FORMATt  3F 1 0 . 4 ) 

FORMAT{  3F7 . 4 ) 

FORMAT( 2F10 .4, 12) 

FORMAT( 13 ) 

FORMAT(F10.5) 

FORMAT( 3F 1 5 . 10) 

FORMAT(2X, 'INPUT  PIVOT  LOCATION:  XP/C=  ') 

FORMAT(F6.4) 

FORMAT{ 2X, 'IS  THERE  A  TRAILING  EDGE  FLAP  [Y/N]?  -  ') 

FORMAT(2X, 'INPUT  POSITIVE  EQUILIBRIUM  YAW  ANGLE  IN  DEGREES:  ') 
FORMAT{2X, 'INPUT  NEGATIVE  EQUILIBRIUM  YAW  ANGLE  IN  DEGREES:  ') 
FORMAT{ /, 2X, ' SET  ANGLE  OF  ATTACK  TO  ', 13, IX, 'DEGREES’ ) 
FORMAT(/,2X, 'PIVOT  LOCATION  -’,F6.4) 

FORMAT( 2X, 'YAW  ANGLE”  ' , 13 , 2X, ’ FORCE- ' ,F8 . 4, 2X, ' TORQUE-' ,F8 . 4 , 
2X, 'Cf^' ,F8.4) 


t 


V  ' 


*  l'  V  »/  • 


FORMAT*  2X,  •  INPUT  TORQUE  DIRECTION  :  '  ,/, 

2X, 'POSITIVE-RESTORING,  NEGATI VE*DE STAB ILI ZING  [P/N]?  -  ') 
FORMAT* 2X, 'TYPE  Y  TO  CONTINUE  -  ') 

FORMAT* 2X, 'END  OF  PROGRAM',/) 

FORMAT ( 2X , ' DO  YOU  WISH  TO  RUN  THE  PROGRAM  AGAIN  [Y/N]?  -' ) 

READ  FILE  CONTAINING  OPERATING  PARAMETERS 

M=0 . 0 
B=0.0 
CALL  HOME 
WRITE* 5, 5000) 

PAUSE 

WRITE* 5 ,5005) 

READ(5,5075)PARAM 
CALL  OPEN ( 6 , PARAM, 2 ) 

READ* 6, 60 00)  XNAME , XDATE 
READ(6,6005)  XCONF 
READ*  6,6010)  B,C,T 
READ* 6, 60 1 1 )XTW, TWT, TWW 
READ(6, 6020 )  THETA, XRAKE,NTAPS 
READ(6,6030)  NPTS 
DO  10  I«1,NTAPS 
READ* 6, 6040 )Y(  I) 

ENDFILE  6 

GET  PARAMETERS  FROM  TERMINAL 

WRITE*  5,5019) 

READ(5,5021 )NRUN 
WRITE* 5, 5020) 

READ*  5,5030)  QI,IS 
WRITE* 5, 5025) 

READ*  5, 5033)  TC 

RHO=. 002522/* 1.0+.00367*TC) 

RMU-. 14248E-8*TC+. 35242E-6 

CALCULATE  QU  -  UPSTREAM  DYNAMIC  PRESSURE 

IF* IS.EQ. 'H' )  GO  TO  30 
QU-1.0898*QI-. 11934 
GO  TO  40 

QU«1.1491*QI-. 73495 
CONTINUE 

VU* SQRT* 2 . 0*QU/RHO ) 

RENU-C* RHO* VU/ ( RMU* 12.) 

CF-C/12. 

BF-B/12. 

TF-T/12. 

AR***.4277 


%  A  *.  %  A  A  A *  % 


DATA  COLLECTION 


WRITE* 5, 5100) 

READ*  5, 50 10 ) IB 

IFflB.NE. 'Y' )  GO  TO  49 
WRITE*  5,5110) 

READ (5,5120)M,BS 
IF(M.EQ.O.O)  GO  TO  48 
WRITE* 5, 7000) 

READ(5,7005)PIVT 
WRITE* 5, 7006) 

READ(5,5010)NT 
WRITE* 5, 7007) 

READ* 5 ,6040) YAP 
WRITE* 5, 7008) 

READ*  5,6040 ) YAN 
DO  75  1-2,22,2 
J-I-2 

WRITE*  5,7010)J 
WRITE* 5, 7004) 

READ*  5, 5010)IC 

IF* IC .NE . ' Y ' )  GO  TO  52 

WRITE(5,7012) 

READ*  5, 50 1 0 )IN 

CALL  ADCONV(NPTS, VBARJ 

F=M*VBAR+BS 

IF*  IN.EQ.  'N*  )  F— F 

TM(I)-F*ARM 

CM* I )*TM( I )/( QU*BF*CF*CF) 

WRITE(5,7018)PIVT 

WRITE*  5, 7020) J , F, TM{  I ) ,CM{ I ) 

CONTINUE 

WRITE  DATA  AND  RESULTS 
L*5 

CALL  HOME 
IP= 'N* 

WRITE* L, 1790) 

WRITE*  L, 1800 )XCONF, XDATE 
WRITE(L, 1801 )XNAME,NRUN 

W  RI TE  ( L ,  1 8  0  5  )  QU ,  Q I ,  RENU ,  TC ,  VU ,  T ,  RHO ,  C ,  RMU 
IF(NT.NE. 'Y' )GO  TO  95 
WRITE*  L,  1806) 

GO  TO  96 
WRITE*  L, 1807) 

CONTINUE 
IF*L.EQ.5)  PAUSE 
WRITE*  L, 1840) 


■  aw  AV|  i  i 


i  a-J  «-> 


WRITE ( L, 1850 )PIVT 
DO  100  I«2,22,2 
J«=I-2 

WRITE(L,  1860  )  J,  TM(  I )  ,CM(  I ) 

WRITE(L, 1870 )YAP, YAN 
IF( L. EQ. 2 )  WRITE{2, 1831 ) 

WRITE ( 5,5072) 

READ( 5 , 50 1 0 ) IP 

IF( IP.NE. * Y' )  GO  TO  9990 

L=2 

GO  TO  91 

FORMAT(5X, ' SYMMETRICAL  FAIRING  SECTIONS' ,//, 5X, 

•2-DIMENSIONAL  TORQUE  COEFFICIENT  TESTING' 

,/, 5X, 'AND  DETERMINATION  OF  HYDRODYNAMIC  CENTER') 

FORMAT(//,5X, 'MODEL  CONFIGURATION  -  • , 40A1 ,/ , 5X , • DATE  -  • 

, 20A1 ) 

FORMAT(//,5X, 'TEST  RUN  NAME  -',6A1,/, 

5X, 'TEST  RUN  NUMBER:  ',13) 

FORMAT(36X, 'TRAILING  EDGE  FLAP:  YES') 

FORMAT(36X, 'TRAILING  EDGE  FLAP:  NONE') 

FORMAT (///,5X, 'POSITIVE  EQUILIBRIUM  YAW  ANGLE  =  *,F7.2, 

'  (DEGREES)',/, 

5X. 'NEGATIVE  EQUILIBRIUM  YAW  ANGLE  ■=  ',F7.2,'  (DEGREES)’) 

FORMAT( 1H1 ) 

FORMAT(//,5X, 'PIVOT' ,5X, 'YAW* ,7X, 'TORQUE' ,4X, 'TORQUE' ,/, 

5X, 'POINT' ,5X, 'ANGLE* , 15X, 'COEFFICIENT' ,/, 

5X, ' (XP/C) ' ,4X, ' (DEG) ' ,5X, ' (FT-LBS) ' ,/, 

5X , 4 1 v ' -  '  )  ) 

FORMAT(4X,F6.4) 

FORMAT ( 16X,I3,5X,F8.4,4X,F8.4) 

FORMATf  /// , 5X , '  QU*= '  ,  F8 . 4 ,  IX,'  (PSF)  '  ,  14X,  'QI='  ,F8.4,  IX,  '  (PSF)  '  ,/, 
5X, 'UPSTREAM  REN=  ', E8 . 3 , 9X, ' TEMPERATURE= ', F6 . 2 , IX ,'( CENTIGRADE ) 
’,/, 

5X,  'VEL*' ,F10.3, IX,  '  (FT/ SEC)  '  ,8X,  '  THICKNESS'  ,  F6. 2,  IX,  ’  ( IN . )  '  ,/ , 
5X, ' RHO*' , FI 0.8, IX, ' (TE  UNITS) ' ,6X, 'CHORD=' ,F6.2, IX, ' (IN. ) ' ,/, 

5X , '  MU= ' , F 1 2 . 10, IX, ' (TE  UNITS) ' ,5X, ' SPAN=’ ,F6.2, IX, ' (IN. ) ' ) 

WRITE( 5, 9001 ) 

READ{ 5, 50 10 ) IA 

IF( IA.EQ. ' Y' )  GO  TO  1 

WRITE( 5 , 9000 ) 

STOP 

END 
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